Draft version February 4, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 2/19/04 



THE INFRARED LUMINOSITY FUNCTION OF GALAXIES AT REDSHIFTS z = l AND z ~ 2 IN THE 

GOODS FIELDS 

K. I. CAPUTI^'^ G. Lagache\ Lin Yan-\ H. Dole\ N. Bavouzet\ E. Le Floc'h*, P. L CHor\ G. Helou^ and N. 

Reddy'' 

Draft version February 4, 2008 

ABSTRACT 

We present the rest-frame 8 /im luminosity function (LF) at redshifts z = \ and ~ 2, computed from 
Spitzer 24 /im -selected galaxies in the GOODS fields over an area of 291 arcmin^. Using classification 
criteria based on X-ray data and IRAC colours, we identify the AGN in our sample. The rest-frame 
8/imLF for star-forming galaxies at redshifts z ~ 1 and ~ 2 have the same shape as at z ~ 0, 
but with a strong positive luminosity evolution. The number density of star-forming galaxies with 
log]^o(^-^E''™) > 11 increases by a factor > 250 from redshift z ~ to 1, and is basically the same at 
z — 1 and ^ 2. The resulting rest-frame 8 /im luminosity densities associated with star formation at 
z — I and 2 are more than four and two times larger than at z 0, respectively. We also compute the 
total rest-frame 8 fim LF for star- forming galaxies and AGN at z ~ 2 and show that AGN dominate 
its bright end, which is well-described by a power-law. Using a new calibration based on Spitzer 
star-forming galaxies at < z < 0.6 and validated at higher redshifts through stacking analysis, we 
compute the bolometric infrared (IR) LF for star- forming galaxies at z = 1 and ^ 2. We find that 
the respective bolometric IR luminosity densities are (1.2 ± 0.2) x 10^ and (6.61} q) x 10* LqMpc^'^, 
in agreement with previous studies within the error bars. At z ~ 2, around 90% of the IR luminosity 
density associated with star formation is produced by luminous and ultraluminous IR galaxies (LIRG 
and ULIRG), with the two populations contributing in roughly similar amounts. Finally, we discuss 
the consistency of our findings with other existing observational results on galaxy evolution. 

Subject headings: galaxies: luminosity function - infrared: galaxies - galaxies: high-redshift - galaxies: 
evolution 



I. INTRODUCTION 

Since the Spitzer Space Telescope (Werner et al. 2004) 
became operational in December 2003, very important 
progress has been made in understanding the nature 
and properties of infrared (IR) galaxies. This progress 
has been revolutionary, in particular, for the study of 
galaxies at high redshifts (z > 1), to which all of the 
previous IR facilities operating in the wavelength range 
A ~ 5 — 200 fim. had basically no access. Previous mis- 
sions as the Infrared Astronomical Satellite (IRAS) and 
the Infrared Space Observatory (ISO) allowed to pro- 
duce multiple studies of mid- and far-IR galaxies, but 
they were restricted to lower redshifts (z < 1) due to 
their sensitivity limits. Until the launch of Spitzer, our 
vision of the high-redshift IR Universe was biased to 
the relatively small number of galaxies detected in sub- 
millimetre and millimetre surveys (e.g. Scott et al. 2002; 
Webb et al. 2003; Greve et al. 2005). 

The sensitivity achieved by the Multiband Imaging 
Photometer for Spitzer (MIPS; Rieke et al. 2004) at 
24 /im is enabling for the first time to conduct system- 
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atic studies of IR galaxies at high redshifts. Several re- 
cent works have shown that, in contrast to what hap- 
pens in the local Universe, the IR extragalactic light is 
increasingly dominated by luminous and ultra-luminous 
IR galaxies (LIRG and ULIRG, respectively) with in- 
creasing redshift (e.g. Le Floc'h et al. 2004; Lonsdale et 
al. 2004; Yan et al. 2004; Le Floc'h et al. 2005; Caputi 
et al. 2006a,b). These LIRG and ULIRG constitute an 
important fraction of the most massive galaxies present 
at z > 1 (Caputi et al. 2006b). 

In a minor but non-negligible fraction of high-redshift 
IR galaxies, the IR emission is produced by the pres- 
ence of an active galactic nucleus (AGN). The exact pro- 
portion of AGN-dominated IR galaxies is actually not 
known, and the determination of such ratio is one of the 
main problems of IR astronomy. A definitive AGN/star- 
forming galaxy separation requires the knowledge of the 
far-IR spectral energy distribution (SED) of these galax- 
ies. Unfortunately, this is not possible for most of high-z 
galaxies, as their far-IR emission is usually below the 
confusion limits at far-IR wavelengths (Dole et al. 2004). 
This separation is also complicated by the existence of 
mixed systems, where both star-formation and AGN ac- 
tivity significantly contribute to the IR emission (e.g. 
Lutz et al. 2005; Yan et al. 2005; Le Floc'h et al. 2006). 
However, the AGN discrimination is essential to disen- 
tangle how much of the IR energy density is associated 
with star-formation. 

The study of a galaxy luminosity function at different 
redshifts allows to understand the composition of the ex- 
tragalactic background as a function of look-back time. 
The analysis of the changes of the LF with redshift is one 
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of the most direct methods to explore the evolution of a 
galaxy population. The first studies of the IR-galaxy LF 
in the local Universe and at low (z < 1) redshifts have 
been based on IRAS and ISO data (e.g. Saunders et al. 
1990; Xu 2000; Takeuchi, Yoshikawa & Ishii 2003; Pozzi 
et al. 2004; Serjeant et al. 2004; Takeuchi et al. 2006). 
Using the most recent Spitzer/MIPS data, Le Floc'h et 
al. (2005) analysed in detail the evolution of the IR LF 
from z = to z ^ 1. They found a positive evolu- 
tion both in luminosity and density between these two 
redshifts, implying that IR galaxies were more numerous 
and the IR output was dominated by brighter galaxies at 
z ~ 1 than at z ~ 0. The IR galaxy LF at higher redshifts 
have been explored by other authors (Pcrcz-Gonzalcz et 
al. 2005; Babbedge et al. 2006). 

Rest-frame 8 /xm luminosities, in particular, are of main 
relevance for star-forming galaxies as they contain in- 
formation on polycyclic aromatic hydrocarbon (PAH) 
emission. PAH molecules characterise star-forming re- 
gions (Desert, Boulanger & Puget 1990) and the associ- 
ated emission lines dominate the SED of star-forming 
galaxies between wavelengths A = 3.3 and 17 /im, 
with a main bump located around 8 /im . Rest-frame 
8 fim luminosities have been confirmed to be good indi- 
cators of knots of star formation (Roussel et al. 2001; 
Forster-Schreiber et al. 2004; Calzetti et al. 2005) and of 
the overall star-formation activity of star-forming galax- 
ies (e.g. Wu et al. 2005), except in low-luminosity galax- 
ies with intense ultraviolet (UV) radiation fields (Gal- 
liano et al. 2005). 

In this work we compute the rest- frame 8 /im LF at 
redshifts z = 1 and z 2, using 24 /im -selected galax- 
ies in the two fields of the Great Observatories Ori- 
gins Deep Survey (GOODS; Giavalisco et al. 2004). At 
z ~ 2, where the fraction of AGN appears to be sig- 
nificant, we analyse separately the LF for star- forming 
galaxies and for the total IR-galaxy population. The 
two GOODS fields cover a smaller area than those anal- 
ysed by some other previous studies of the IR LF. How- 
ever, they benefit from uniquely deep homogeneous pho- 
tometric datasets, ranging from the X-rays to radio wave- 
lengths, as well as an important spectroscopic coverage. 
As we explain in Section [21 this makes possible an al- 
most complete identification of 24 fj,m galaxies down to 
faint fluxes and the derivation of accurate redshift deter- 
minations (cf. also Caputi et al. 2006a, c). These two 
characteristics are essential for a proper computation of 
the LF at high redshifts, without any conclusion relying 
either on completeness or selection-function corrections. 

The layout of this paper is as follows: in Section [2l 
we describe in detail the selection of our 24 /im -galaxy 
samples at redshifts 0.9 < z < 1.1 and 1.7 < z < 2.3. 
In Section [3l we explain how we perform the separation 
between star-forming galaxies and AGN within our sam- 
ple. We compute the rest-frame 8/imLF at z = 1 in 
Section |4] and analyse its evolution from z 0. In Sec- 
tion \5\ we present the rest-frame 8 /^m LF at z ~ 2 and 
extend the analysis of the evolution up to this high red- 
shift. Later, in SectionlHl we use a new empirical calibra- 
tion based on Spitzer galaxies to obtain the bolometric 
IR LF at different redshifts. Finally, in Sections [7] and 
[51 respectively, we discuss our results and present some 
concluding remarks. We adopt throughout a cosmology 



with Ho = 70kms-iMpc~\ = 0.3 and 17a = 0.7. 

2. THE IR GALAXY SAMPLE IN THE GOODS FIELDS 

The GOODS fields, namely the GOODS/Chandra 
Deep Field South (GOODS/CDFS) and 
GOODS/Hubble Deep Field North (GOODS/HDFN), 
have been observed by Spitzer as one of the cycle- 1 
Legacy Science Programs (P.I. Mark Dickinson). Ex- 
tended areas of the CDFS and HDFN have also been 
observed as part of the Spitzer IRAC and MIPS Guar- 
anteed Time Observers (GTO) programs (P.I. Giovanni 
Fazio and George Rieke, respectively). 

GOODS/IRAC maps at 3.6 to 8/imand MIPS maps 
at 24/imare now publicly available. The corresponding 
GOODS pubhc 24 /xm catalogues have been constructed 
using prior positional information from the IRAC 3.6 and 
4.5 /im images and by an additional blind extraction of 
24 /im sources. The resulting 24 /im catalogues are basi- 
cally reliable and complete for galaxies with fluxes down 
to 5(24 /xm) = 80 /iJy (Chary et al., in preparation, and 
cf. the Spitzer GOODS website®). For a comparison, 
we note that the 24 /im catalogue constructed from the 
shallower MIPS/GTO observations of the CDFS achieves 
^ 80% completeness and only has ~ 2% of spurious 
sources at a similar flux level (Papovich et al. 2004). Al- 
though, in principle, fainter sources can be detected in 
the deeper GOODS images, we decide to only use the 
conservative GOODS 5(24 /im) > 80/iJy galaxy cata- 
logues for the selection of our 24 /im -galaxy samples at 
z ~ 1 and z ~ 2. In this way, our computed LF are 
virtually not affected by incompleteness corrections (cf. 
sections [H and |5|) . 

2.1. Multiwavelength analysis and redshift 
determinations for 24 /im sources in the 
GOODS/CDFS 

In the GOODS/CDFS, we restrict our analysis to the 
131 arcmin^ which have deep J and i^s-band cover- 
age by the Infrared Spectrometer and Array Camera 
(ISAAC) on the 'Antu' Very Large Telescope (Antu- 
VLT) (GOODS/EIS vl.O release; Vandame et al., in 
preparation). We used the Kg < 21.5 (Vega mag) 
galaxy catalogue constructed by Caputi et al. (2006c) 
to identify the 24 //m galaxies in the GOODS/CDFS cat- 
alogue, using a matching radius of 2". The percent- 
age of 24 /im galaxies with double iiTg-band identifications 
within this radius is only < 8% and 95% of the as- 
sociations can be done restricting the matching radius 
to 1.5" (Caputi et al. 2006b). In all cases of multi- 
ple identifications, we considered that the counterpart 
to the 24 /im source was the Kg galaxy closest to the 
24 /im source centroid. The Kg < 21.5 mag catalogue 
allows us to identify 515 24 /im galaxies within the 131 
arcmin^ area, i.e. ~ 94% of the 24 /xm galaxies with 
5(24 /im) > 80/iJy in this field. 

Caputi et al. (2006c) measured multiwavelength pho- 
tometry for all their Kg < 21.5 mag galaxies. They ran 
SEXTRACTOR (Bertin & Arnouts 1996) in 'double- 
image mode' to perform aperture photometry on the 
GOODS/EIS vl.O J-band images, centred at the posi- 
tion of the i^s-band extracted sources. They also looked 

^ http:/ /data. spitzer. caltech.edu/popular/goods 
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for counterparts of the Ks < 21.5 mag sources in the pub- 
lic GOODS Advanced Camera for Surveys (ACS) cata- 
logues, which provided photometry in the B, V, 1775 and z 
bands. The stellarity parameter measured on the 2;-band 
images allowed them to separate out galactic stars. Fi- 
nally, they ran SEXTRACTOR on the Spitzer/IRAC 3.6 
and 4.5 /im images to identify the Kg < 21.5 mag galaxies 
and measured aperture photometry at these longer wave- 
lengths. We refer the reader to Caputi et al. (2006c) for 
additional details about the photometric measurements 
and applied aperture corrections. 

Caputi et al. (2006c) obtained an estimated redshift 
for each one of their galaxies modelling their stellar 
SED from the B through the 4.5 /im bands. They used 
the pubhc code HYPERZ (Bolzonella, Miralles & Pello, 
2000) with the GISSEL98 template library (Bruzual & 
Chariot 1993) and the Calzetti et al. (2000) reddening 
law to account for internal dust extinction. 

The HYPERZ redshift estimates have been replaced 
by C0MB017 photometric redshifts (Wolf et al. 2004) 
for those galaxies with magnitudes R < 23.5 mag at 
redshift z < 1, which is the regime of higher accuracy 
for C0MB017. In these cases, the SED fitting has 
been constrained to the C0MB017 redshifts. The cross- 
correlation of the GOODS/CDFS 24 /y,m catalogue with 
the Caputi et al. (2006c) Kg < 21.5 mag catalogue di- 
rectly gives us estimated redshifts and best-fitting SED 
models for all the identified 24 /im galaxies. 

2.2. Multiwavelength analysis and redshift 
determinations for 24 /xm sources in the 
GOODS/HDFN 

In the GOODS/HDFN, we followed a similar strat- 
egy for the analysis of sources as in the GOODS/CDFS. 
However, unfortunately, we only have access to deep Kg- 
band data for a part of this field (Reddy et al. 2006b). 
Thus, we used the Spitzer/IRAC S.Qi^m maps to iden- 
tify the 24 /xm galaxies. We analysed in this case the 
entire GOODS/HDFN region, i.e. the 160 arcmin^ with 
deep GOODS/HST-ACS coverage. We ran SEXTRAC- 
TOR on the IRAC 3.6 and 4.5 /im images. We con- 
structed a catalogue of 3.6 /im sources, accepting only 
those objects also identified in the 4.5 /im band. To en- 
compass the technique applied by Caputi et al. (2006c) 
on the IRAC maps of the GOODS/CDFS, we measured 
photometry in circular apertures of 2.83"-diameter'' and 
applied aperture corrections of 0.50 and 0.55 mag to 
the 3.6 and 4.5 /tm magnitudes, respectively. We then 
used this 3.6 /xm catalogue to identify the 24 /xm sources 
in the GOODS/HDFN, using a matching radius of 2". 
This allows us to identify 856 24 /im galaxies in the 160 
arcmin^ of the GOODS/HDFN, i.e. - 95% of the 
24 /xm galaxies with 5(24 /xm) > 80 /tJy in this field. The 
identification completeness achieved for 24 /xm galaxies in 
this field using 3.6 /xm sources is similar to the identifi- 
cation completeness obtained for 24 /xm galaxies in the 
GOODS/CDFS using if^-band sources. This indicates 
that the two identification methods are basically equiv- 
alent. In any case, the IRAC 3.6 and 4.5 /tm data are 
incorporated in the SED modelling of all the sources in 
the two fields. 

The aperture size has been chosen in correspondence to the 
aperture sizes used in the GOODS ACS catalogues. 



We followed up in the optical bands those IRAC 
3.6 /xm objects which were counterparts to 5(24 /xm) > 
80 /tJy sources. Once more, we used the public GOODS 
ACS catalogues to obtain aperture photometry in the 
B, V, Irr5 and z bands. In addition, we looked for coun- 
terparts of the 3.6 /xm sources in the U and HK'-band 
images of the GOODS/HDFN (Capak et al. 2004). Al- 
though these images are relatively shallower than the 
other optical/near-IR data available for this field, we de- 
cided to include these data to improve the SED coverage. 
Finally, we incorporated the deep J and Kg-hand data 
from Reddy et al. (2006b) for those galaxies lying in the 
region where these data were available (< 40% of the 
analysed area). 

We used the multiwavelength data from the U to 
the 4.5 /tm bands to model the SED and obtain pho- 
tometric redshifts for all of our 24 /xm galaxies in the 
GOODS/HDFN using HYPERZ, in an analogous way to 
that in Caputi et al. (2006c). As in the latter, we applied 
a set of criteria to control the HYPERZ output: 1) the 
photometric redshifts for galaxies detected in the shal- 
low U-band catalogues were constrained to a maximum 
value Zphoi = 2, as bright U-band sources are unlikely to 
be beyond these redshifts; 2) analagously, the estimated 
redshifts of galaxies not-detected in the U-band but de- 
tected in the B-band were constrained to a maximum 
value Zphot = 4; 3) for the GOODS/HDFN catalogue, we 
found that HYPERZ produced an ovcrdcnsity of galax- 
ies in the redshift range 1.5 < ZphoL < 1-7. Compar- 
ison with spectroscopic redshifts (see below) suggested 
that this overdensity was an artefact of HYPERZ ap- 
plied to our sample. Thus, to test these possible spurious 
redshifts, we double-checked the fitting of all the galax- 
ies with HYPERZ redshift 1.5 < Zphot < 1-7 using the 
PEGASE library (Le Borgne & Rocca-Volmerange 2002). 
We kept the HYPERZ solution for those sources con- 
firmed by PEGASE as belonging to the 1.5-1.7 redshift 
range. For all the remaining 1.5 < Zphot < 1-7 galaxies, 
we replaced the photometric redshift by the PEGASE 
estimate. This strategy improved the agreement with 
spectroscopic redshifts. The percentage of galaxies with 
PEGASE redshifts in our final 24 /xm catalogue for the 
GOODS/HDFN is 5%. 

2.3. The final IR galaxy samples in the combined 
GOODS fields 

Our final 24 /xm catalogue contains 1371 24 /xm sources 
with S'(24 /tm) > 80 /iJy over a total area of 291 arcmin^. 
We identified only 22 out of 1371 sources as galactic stars. 
All the remaining sources are galaxies. Our aim is to 
separate two sub-samples of galaxies from this final cat- 
alogue: 1) the 24 /tm galaxies with redshifts 0.9 < ;z < 1.1 
for the computation of the IR LF at z = 1; 2) the 
24 /xm galaxies with redshifts 1.7 < 2; < 2.3 for the com- 
putation of the IR LF at z ~ 2. 

We performed a final step before separating the two 
definitive sub-samples of 24 /xm galaxies used in this 
work. In addition to the wealth of photometric data, 
both GOODS fields benefit from an important amount 
of spectroscopic data, most of which are publicly avail- 
able (Cohen et al. 1996; Le Fcvre et al. 2004; Wirth 
et al. 2004; VanzeUa et al. 2005, 2006; Choi et al., in 
preparation; among others). Some additional redshifts 
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in the GOODS/CDFS have been kindly made available 
to us by Frangois Hammer and Hector Flores. We com- 
piled these data and found that more than 45% of our 
24 /im galaxies in the combined fields had spectroscopic 
rcdshifts. We incorporated these spectroscopic redshifts 
into our catalogue, which superseded the corresponding 
photometric values. The finally discarded photometric 
redshifts have been used to assess the quality of our red- 
shift estimates. Figure [T] shows the comparison between 
photometric and spectroscopic redshifts for the galax- 
ies in our sample for which both redshifts are available. 
We observe a good agreement between photometric es- 
timates and real redshifts. The distribution of relative 
errors dz — [zphot ~ 2spec)/(l + Zspec) has a median value 
-0.007 and the dispersion is az = 0.05. 

From the definitive redshift catalogue which in- 
corporates spectroscopic redshifts, we select those 
24 ^m galaxies lying at 0.9 < z < 1.1 and 1.7 < z < 2.3. 

The 0.9 < z < 1.1 sample is composed of 227 galaxies 
with 5(24 /im) > 80/iJy and a median redshift z = 1.00. 
We use this sample to compute the IR LF at z — 1 . More 
than 60% of these galaxies have spectroscopic redshifts 
2spcc- The quality of photometric redshifts is similar as 
that for the total sample: the median of relative errors 
dz = {zphot — Zspec) /(I + Zspec) IS -0.01 and the dispersion 
is (Tz = 0.05. In the computation of the IR LF at z = 1, 
we consider that these errors only affect those galaxies 
with photometric redshifts (< 40%). 

Our 1.7 < z < 2.3 sample contains 161 24 /xm galaxies 
with S'(24 /im) > 80 /iJy. This is the sample we use to 
compute the IR LF at redshift z 2. The median red- 
shifts of these 161 galaxies is z — 1.93. Although for 
practicity we refer to these galaxies as the z ~ 2 sample, 
all the calculations made in Sections [5] and [6] take into 
account the actual median redshift value. More than 
15% of the galaxies selected with 1.7 < z < 2.3 have 
spectroscopic redshifts. The quality of photometric red- 
shifts for the z 2 sample can also be assessed from 
figure [TJ We see that the agreement between photomet- 
ric and spectroscopic redshifts is still very reasonable for 
this high-redshift sample. The distribution of relative 
errors dz = (zphot - 2spec)/(l + Zspec) has a median - 
0.01 and a dispersion az = 0.06. This statistics has been 
computed based on all sources (i.e. AGN included; see 
below). This suggests that the SED templates we use 
to derive photometric redshifts are suitable for all our 
sample. The photometric redshift error bars affect the 
majority of galaxies in our z ^ 2 sample and are taken 
into account in the computation of the corresponding LF, 
as we explain in Section [5l 

We note that the galaxies with spectroscopic redshifts 
are representative of our entire 24 /im sample in each 
of the considered redshift bins (0.9 < z < 1.1 and 
1.7 < z < 2.3). The two panels in Figure [2] show the 
rest- frame 8 /im luminosities of all of our galaxies (empty 
histograms) and those of galaxies with spectroscopic red- 
shifts (shaded histograms), at these different redshifts. 
Details on the calculation of 8 /im luminosities are given 
in Section 31 From Figure [H we can see that galax- 
ies with spectroscopic redshifts basically span the whole 
range of IR luminosities considered in this work. Thus, 
the errors derived from the comparison of photometric 
and spectroscopic redshifts are applicable to the entire 
IR LF. 



3. THE NORMAL/ACTIVE GALAXY SEPARATION 

In this work, we would like to compare the IR LF for 
star- forming galaxies only with the total IR LF. To do 
this, we need to identify the active galaxies present in 
our sample. 

One of the most efficient ways of identifying AGN is 
through their X-ray emission. The GOODS fields have 
deep X-ray coverage obtained with the Chandra X-ray 
Observatory: the IMs maps for the CDFS (Giacconi et 
al. 2002) and the 2Ms maps for the HDFN (Alexander et 
al. 2003). We used the corresponding public X-ray cata- 
logues to identify the AGN within our sample. However, 
given the depth of these catalogues (especially that of the 
HDFN) , X-ray sources include not only quasars and AGN 
but also powerful starburts which also emit in X-rays. To 
separate the two classes of X-ray sources, an optical ver- 
sus X-ray fiux diagram can be used. Figure [3] shows the 
i?-band magnitude versus the soft X-ray flux of the X- 
ray-detected galaxies in our 24 /im sample in the HDFN. 
The i?-band magnitudes of our galaxies have been inter- 
polated using the V and 1775 magnitudes. This plot is 
similar to that presented in Alexander et al. (2003). The 
dashed line shows the empirical separation between nor- 
mal galaxies and AGN, as calibrated by Hornschemeier 
et al. (2001). Using this diagram, we identify the X-ray 
detected AGN within our 24 /im sample and, in particu- 
lar, those at 0.9 < z < 1.1 and 1.7 < z < 2.3. 

Some AGN with weak soft X-ray fluxes but signiflcant 
emission in the hard bands can contaminate the normal 
galaxy region in the i?-band magnitude versus soft X- 
ray flux diagram. These AGN are characterised by a flat 
photon index F < 1.0 (e.g. Hornschemeier et al. 2003). 
We also looked for these kinds of objects to identify the 
AGN present in our sample. 

It is known, however, that the X-ray selection can be 
incomplete for the selection of AGN. Other active galax- 
ies exist, which are not detected even in deep X-rays sur- 
veys. A complementary method to select active galaxies 
can be developed based on the analysis of the IR colour 
excess in the Spitzer IRAC bands. Figure 2] shows the 
(3.6-8 /im) versus (5.8-8 /im) colours for all the galaxies 
with redshift z > 1.5 in our 24 /im sample. Empty circles 
correspond to all those galaxies not classified as AGN 
using X-ray data (either not detected in X-rays or X-ray 
sources classified as starbursts). Filled squares indicate 
the X-ray classified AGN. We restrict this diagram to 
high redshift sources for the following reason. The stel- 
lar bump centred at rest- frame wavelength A ~ 1.6 /xm 
is shifted into the IRAC bands at z > 1.5. For active 
galaxies, the galaxy SED at the same rest-frame wave- 
lengths is dominated by a power-law continuum. Thus, 
it is expected that an IRAC-based colour-colour diagram 
is able to separate the AGN through their IR excess. At 
low redshifts, this separation is much less clear, espe- 
cially because star-forming galaxies with PAH emission 
can mimic the IR excess. Similar colour-colour plots have 
been used with the purpose of separating normal and ac- 
tive galaxies elsewhere (e.g. Lacy et al. 2004; Stern et 
al. 2005; Caputi et al. 2006b). 

Inspection of figure |4] shows that X-ray-selected AGN 
display a wide range of (3.6-8 /im) and (5.8-8 /im) colours, 
while the vast majority of "normal galaxies" (i.e. non 
X-ray-classified AGN) appear on the left-hand side of 
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this diagram, with a colour (5.8-8 ^m) < 0.2 (AB). As 
we mentioned above, the relatively blue colours are pro- 
duced by the stellar SED bump mapped at the IRAC 
wavelengths. The galaxies lying on the right-hand side, 
on the contrary, present an excess in the SED contin- 
uum which is characteristic of AGN. Thus, based on this 
diagram, we adopt an empirical colour cut to produce 
an additional AGN selection criterion: all the z > 1.5 
galaxies with (5.8-8 /im) > 0.2 (AB) within our sample 
are classified as AGN. This same additional AGN selec- 
tion criterion has been used by Caputi et al. (2006b). 

We would like to note that, while this colour cut pro- 
duces a safe criterion to select additional active galaxies, 
it is possibly not complete. The dispersion of colours dis- 
played by X-ray-selected AGN suggests that other active 
sources -not detected in X-rays and with no IRAC colour 
excess- could also exist among the 24 /im galaxies. On 
the other hand, some of the X-ray classified AGN could 
be composite systems, where a fraction of the bolometric 
IR luminosity is actually due to star formation. Unfor- 
tunately, no AGN selection criterion appears to be both 
complete and reliable at the same time (cf. e.g. Barmby 
et al. 2006). As we do not have information on the far- 
IR emission of our galaxies, our separation criteria are 
possibly the most adequate to discriminate AGN. 

For our sample of 227 24 /im galaxies with redshift 
0.9 < z < 1.1, only the X-ray criteria have been applied. 
We identify 23 out of 227 galaxies as AGN, i.e. ~ 10% of 
the sample. We will exclude the AGN from our sample 
in order to determine the IR LF for star-forming galaxies 
Bit z = 1, but we note that the inclusion of AGN only has 
a minor impact on the LF at this redshift. 

For the sample at redshifts 1.7 < z < 2.3, we ap- 
plied both selection criteria to separate AGN (X-ray 
and IRAC-colour classifications). The fraction of active 
galaxies at these redshifts appears to be more impor- 
tant than at z 1. We identify 29 AGN among our 161 
24 /im galaxies at 1.7 < z < 2.3, i.e. ~ 18% of the sample. 
23 out of these 29 AGN have been identified using X-rays 
and the remaining 6 AGN have been classified through 
their IRAC colours. As we will see below, the LF com- 
puted including and excluding AGN have non-negligible 
differences, because these objects dominate the bright 
end of the IR LF at these high redshifts. Throughout 
this paper, when we refer to the star-forming galaxies 
at redshift z ~ 2, we mean our sample of 161-29=132 
objects which we have not classified as AGN at these 
redshifts. 

4. THE REST-FRAME 8 fimhF AT REDSHIFT z = l 
4.1. The k-corrections from 11.4 — 12.7 to 8 fiia 

Before computing the rest-frame 8 /xmLF at z ~ 2, we 
aim to understand its evolution from z ~ to redshift 
z = 1. For this, we compute the rest-frame 8/imLF for 
our 204 24 /im -selected star- forming galaxies in the red- 
shift range 0.9 < z < 1.1. AGN have been excluded from 
this analysis. AGN constitute ~ 10% of our sample with 
0.9 < z < 1.1 and their exclusion does not significantly 
change the shape of the 8 iim LF at z = 1 . This is in con- 
trast to what we find at z ~ 2, where AGN constitute a 
somewhat higher fraction of sources which dominate the 
bright end of the rest-frame 8 /xm LF (cf. Section [5|) . 

We compute the rest-frame 8 /xm luminosity (i^i^^'") 
of each galaxy as 



i/L^pm ^ i.47rfc(Arf)S'(24/im)d£(z), where 5(24^m) is 
the 24 /im flux, d-!^{z) is the luminosity distance and 
fc(Ai.f), the corresponding k-correction at the rest-frame 
wavelength Aif. The width of the redshift bin we con- 
sider 0.9 < z < 1.1 implies that the observed 24 /im maps 
rest-frame wavelengths 11.4 < Ajf < 12.7/im. We need 
then to apply k-corrections to convert the rest-frame 
11.4— 12.7/im into 8/j,mfluxes. 

To compute these k-corrections, we analyse differ- 
ent sets of IR galaxy templates available in the liter- 
ature, namely the models by Chary & Elbaz (2001) 
and Elbaz et al. (2002); Dale et al. (2001) and Dale 
& Helou (2002); and Lagache et al. (2004). We con- 
volve the SED templates in all these models with the 
transmission function of the 24 /im filter and obtain the 
relation between the fluxes at 11.4 — 12.7 and 8 /im. Fig- 
ure[5]shows the A-to-8 /im k-corrections in the wavelength 
range Art = 11.4— 12.7 /zm. Different line styles indicate 
the k-corrections obtained with different SED templates. 
The solid and dotted lines correspond to the range of 
k-corrections derived for galaxies with bolometric IR lu- 
minosities LiR > 10^^ Lq, from the Lagache et al. and 
Chary & Elbaz models, respectively (with thick lines in- 
dicating the median values). The dashed lines show the 
k-corrections obtained with the Dale et al. SED model 
with parameters a = 1.1 and 1.4 (cf. Dale et al. 2001). It 
is clear from inspection of figure [S] that the k-corrections 
between 11.4 — 12.7 and 8 /im obtained with these differ- 
ent models have some significant dispersion. These dif- 
ferences are produced by the limited knowledge on PAH 
emission when modelling the PAH-dominated region of 
a star-forming galaxy SED. 

In this work, we adopt the median k-corrections ob- 
tained with the Lagache et al. (2004) models of star- 
forming galaxies with bolometric IR luminosities Lir > 
10^^ L0 (thick solid line in [5]). As we show in Section 
16.1.21 the Lagache et al. templates produce an 8 /zm - 
to-bolometric IR luminosity conversion quite close to 
that measured on the observed SED of Spitzer galaxies 
(Bavouzet et al. 2006). This suggests that these tem- 
plates incorporate an adequate modelling of the PAH 
emission region in the star- forming galaxy SED. 

4.2. The l/Vmax method 

We compute the rest-frame 8 /xm LF for star- forming 
galaxies at redshift z = 1 using the 1/Vmax method 
(Schmidt 1968). For this calculation, we consider the 
204 star-forming galaxies with redshift 0.9 < z < 1.1 
within our sample. The advantage of the 1/VJnax tech- 
nique is that it allows to compute the LF directly from 
the data, with no parameter dependence or model as- 
sumption. Besides, the normalization of the LF is di- 
rectly obtained from the same calculation. The comov- 
ing volume V^nax — K„ax ~ — 0-9) ^^r each source 
corresponds to the maximum redshift Zmax at which it 
would be included in the catalogue, given the limiting 
flux S'(24 /im) = 80 /iJy, and provided that this redshift 
is smaller than the maximum of the considered redshift 
bin (in this case z ~ 1.1). Otherwise, T^nax is equal 
to the volume corresponding to the 0.9 < z < 1.1 bin 

^^ax — ^^in • 

As we explained in Section [2 the GOODS 
24 /im catalogues are basically complete down to the lim- 
iting flux and, thus, no sample completeness corrections 
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are needed for our catafogues. However, we do apply 
completeness corrections to account for the percentage 
(5-6%) of unidentified 24 sources (cf. Section [5]). 
These identification completeness corrections are very 
small and none of the conclusions presented here depend 
on the application of such corrections. 

We present the results of our rest-frame 8 iim LF at 
redshift z = 1 computed with the l/Knax method in Fig- 
ure [S] (upward-pointing triangles). This LF, as well as 
all the other presented in this work, have been computed 
jointly on the GOODS/CDFS and GOODS/HDFN. Al- 
though we have checked the consistency within the er- 
ror bars of the LF obtained in the two fields sepa- 
rately, the sample variance effects are more important 
than when considering both fields combined (cf. Fig- 
ure [1]). We show the 8/xmLF function computed with 
the 1/Vmax method only in the completeness region of 
8 /zm luminosities (z/L^'^™ ^ 3x 10^° L©), imposed by the 
flux limits of the 24 ^m survey and the considered red- 
shifts. The total comoving volume probed by our survey 
is 1.3 X lO^Mpc^. 

The error bars for this LF values depend not only on 
the number of sources (Poisson statistics), but also on 
the errors in the photometric redshifts and in the k- 
corrections applied. The errors in the photometric red- 
shifts affect only < 40% of our galaxies at 0.9 < z < 1.1, 
given the high percentage of available spectroscopic red- 
shifts. To account for the errors in the photometric red- 
shifts, we made Monte Carlo simulations of our vLf,'^™ 
catalogues. We produced 1000 simulated catalogues, 
each one with the same number of sources as our orig- 
inal 0.9 < z < 1.1 catalogue of star-forming galaxies 
(i.e. 204 sources each) . The redshift of each source has 
been allowed to randomly vary following a Gaussian dis- 
tribution centred at Zccntrc = z — 0.01 and with a dis- 
persion (T^ — 0.05(1 -I- z) (cf. Section [2T3|) . where z is 
the redshift of the source in the original catalogue. The 
redshift of those sources with spectroscopic determina- 
tions have been left unchanged. For the k-corrections, 
we fixed the error bars to e = 0.50, which is roughly 
the dispersion between the different Lagache et al. and 
Dale & Helou model predictions (see Figure [5]) . To in- 
clude these errors in the simulations, we computed the 
rest- frame 8 /^m luminosity lyLf,^"^ of each galaxy in the 
mock catalogue allowing the corresponding k-correction 
to have a random value within the range of its error bar. 
Finally, the LF has been recomputed with the 1/Vniax 
method for each of the mock catalogues, with exactly 
the same procedure as for the original catalogue. From 
the distribution of the LF values in each j/L^'"" bin, we 
determined the error bars on our original l/Knax results. 

For a comparison, we also show the 8/xmLF of star- 
forming galaxies at redshift z ~ (strictly < z < 0.3; 
with median z w 0.2), computed by Huang et al. (2006), 
using the l/Vmax method applied to IRAC 8^mGT0 
data (cross-like symbols in Figure [6]). No error bars 
have been plotted for this LF, as they are significantly 
smaller than the error bars of the LF we determine here. 
The comparison of this 1/Vmax LF with our own de- 
termination at z = 1 shows a substantial increment 
of the density of star-forming galaxies with rest-frame 
8 /xm luminosities log]^Q(z^i^'"") > 10.5, with increasing 
redshift. We note that this behaviour is evident from 
the l/Knax calculation, independently of the parametric 



analysis we discuss below. 

4.3. The maximum likelihood analysis 

The shape of the z ~ LF can be fitted with a 
double-exponential function (Saunders et al. 1990; Pozzi 
et al. 2004; Le Floc'h et al. 2005): 



$(L)dlogi„(L) 



X exp 



dlogio(L), (1) 



where, in this case, L = j/L^'^™. The parameters a and 
1/cr^ correspond to the slopes at the faint and the bright 
ends, respectively. L* is the characteristic I'L^^'^™ lumi- 
nosity where the transition between the faint and bright 
regimes occurs, and $* is the normalization factor. Usu- 
ally, the parameter a is fixed a priori, as the faint-end 
of the LF is poorly constrained. Fixing a = 1.2 (e.g. 
Zheng et al. 2006) and using a minimization tech- 
nique, we obtain that the best-fitting parameters for the 
LF at z ~ are cr = 0.36±0.01, L* = (5.8±0.2) x 10^ L© 
and $* = (5.7 ±0.1) x 10"^ Mpc^^ dex^i. The result- 
ing curve (dotted line in Figure [6|) produces an excellent 
fitting of the l/^max LF at z ~ 0. 

Assuming that the form given in eq.([T]) is also suit- 
able to describe the IR LF for star-forming galaxies at 
higher redshifts, we obtain a second independent calcu- 
lation of the rest-frame 8 /im LF at redshift z = 1 using 
the STY (Sandage, Tammann & Yahil, 1979) maximum 
likelihood (ML) analysis. This is a parametric technique 
which assumes a given shape for the LF. No data binning 
is involved in the calculation. The advantage of the ML 
analysis over the 1/Vmax technique is that the former 
does not contain any assumption on a uniform spatial 
distribution of galaxies. The corresponding likelihood 
estimator reads: 



C[sk\{zi,L^ 



,N 



N 

n 

4=1 



$(sfe,L) dlogio(L) 



(.2) 

where the product is made over the i — 1, N galaxies 
of the sample, ^{s^, L) is the adopted form for the LF as 
a function of the luminosity i, and which depends on the 
parameters su- L\ is the minimum luminosity at which 
the i-th galaxy would be observable, given its redshift Zi 
and the flux limit of the survey. The weighting factors 
Wi allow to take into account completeness corrections 
(Zucca, Pozzetti & Zamorani 1994; Ilbert et al. 2005). 
By maximizing C (or, for simplicity, its logarithm), one 
can obtain the values of the parameters Sk yielding the 
maximum likelihood. The normalization factor $* is re- 
covered after the maximization, by integrating the ob- 
tained maximum likelihood LF without normalization in 
the range of luminosities of the survey, and making it 
equal to the number density of observed galaxies. We 
note that the ML analysis provides a direct calculation 
of the LF (i.e. it does not constitute a fitting procedure 
as the minimization) and is completely independent 
of the LF obtained with the 1/Vmax technique. 

For the case of our rest-frame 8/imLF at z = 1, we 
apply the STY method using eq. H]) and fixing the slopes 
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at the faint and bright ends to the same values as at 
z 0, i.e. a — 1.2 and a — 0.36, respectively. In this 
case, we obtain that the value of the characteristic lumi- 
nosity which maximizes the likelihood estimator is L* = 
j,^*8A^m ^ (3.55+0-5g) X IQIOLq and the derived normal- 



ization factor is = (3.95;[!jg) x IQ-^ Mpc'^ dex"!. 
The error bars on L* have been computed considering 
A(ln£) = —0.5 and the uncertainties derived from the 
Monte Carlo simulations. The degeneracies in parame- 
ter space given by A(ln£) = —0.5 dominate the L* error 
budget. The error bars on have been derived using 
the extreme values of L* (i.e L*± its error). The re- 
sulting curve for the ML LF at z = 1, obtained with a 
double-exponential law with a = 0.36, is indicated with 
a solid line in the upper panel of Figure [HI 

Another possibility is to allow the slope at the bright 
end (l/cr^) to be a free parameter in the ML analy- 
sis. Doing so, we obtain that the ML is produced for: 



IS pr 

(LlOiglD X 10" Lo 



= (2.54^ 

(dotted-dashed line in the upper 



^+0.60^ 
-0.35-' 



-•-0.07' 

and the derived normalization is 
lO-^Mpc-^dex-i 

panel of Figure [6]) . The degeneracy in (cr, L* ) space is 
shown in the lower panel of this same figure. 

The LF obtained with the ML analysis, both in the 
case of a free cr value and fixed a = 0.36, are in 
good agreement with the LF computed with the 1/Vniax 
method. This confirms that the double-exponential law 
in eq. ([T]) also provides a good description of the 8 /im LF 
at high redshifts. The degeneracy existing in the a value 
is due to the limited constraint that our data can put on 
the bright end of the LF at z = 1. In Figure [6l we see 
that the double-exponential forms with a — 0.20 and 0.36 
only differ significantly at the very bright end of the LF 
(i/L^''™ > IO^-'^'^Lq at z = 1). Large-area surveys with 
a significant number of very bright IR galaxies, as for 
example the ^ 2 deg^ Spiizer-COSMOS survey (Sanders 
et al. 2006), are necessary to set tighter constraints in 
(cr, L*) space. 

Finally, we explore whether other functional forms 
could also be suitable to describe the rest-frame 8 /zm LF 
at z = 1. We repeat the calculation of the LF with 
the STY method, but this time using a Schechter (1976) 
function: 



^{L)d\og^^{L) = ^* 



X exp 



L 

L* 



dlogio(i)- 
(3) 

By fixing a — 1.2, we find that the maximum likelihood is 
obtained for a characteristic luminosity L* = lyLl ^ ^™ = 
(7.2j;Q7) X 10^" L0 and the derived normalization is $* — 
(3.88l^;4^) X 10-3 Mpc-3 dex-^. The resuffing ML curve 
is shown with a dashed line in Figure [6] The Schechter 
form actually produces a LF quite close to that obtained 
with the cr — 0.20 double-exponential form, within the 
observed luminosity range of our survey. 

The degeneracy existing in the shape of the IR LF, as 
constrained from our data, produces some uncertainty in 
the determination of the number density of the most lu- 
minous IR galaxies (cf. TableE]). However, as we discuss 
below, this degeneracy has little impact on the derived 
luminosity density. This value is mainly governed by the 
turnover of the LF, which we can properly determine 



here, given the depth of our survey. 

4.4. The evolution of the rest-frame 8 /im LF from z 

io z = 1 

When using the same law to describe the rest-frame 
8/imLF both at redshifts z ~ and z = 1, the differ- 
ences found in the characteristic luminosity L* and the 
normalization parameter $* can be understood as a lu- 
minosity and density evolution: 



L*{Z2 - 1) -i*(zi -0) X 



$*(Z2 = 1) = $*(zi - 0) X 



1 + Z2 
1 + Zl 
1 + Z2 
1 + Zi 



(4) 



where we strictly use zi = 0.2 (the median redshift of the 
Huang et al. sample). 7l and 75 describe the evolution 
of the L* and <&* parameters with redshift. The values 
of these parameters at z ~ and z = 1, corresponding 
in both cases to a double-exponential with a = 0.36, 
produce (cf. Section [ 



7L= 3.5 ±0.4 

75 = -0.7 ±0.1. (5) 

This implies a strong positive luminosity evolution and a 
mild negative density evolution between z and z = 1. 
The mild negative density evolution to z = 1 refers to the 
overall normalization However, it is clear from Fig- 
ure [6] that, within the 8 /xm luminosity range spanned by 
our sample, the density of galaxies at z = 1 is larger than 
that at z ~ 0. This is consistent to what has been found 
by Le Floc'h et al. (2005) from the analysis of the rest- 
frame 15 fim LF and confirms, once more, the increasing 
importance of IR galaxies up to redshift z ^ 1. For 
clarity, the density of galaxies we obtain by integrating 
the rest- frame 8 /im LF above different luminosity cuts at 
different redshifts are shown in Table [5] 

By integrating the LF weighted by the luminosity val- 
ues, over all luminosities, we can obtain the total rest- 
frame 8 /im luminosity density. In fact, for the obtention 
of the total luminosity density, the precise limits of in- 
tegration are irrelevant, provided the turnover of the LF 
is completely contained within these limits. Moreover, 
the use of any of the different laws which are suitable to 
describe the LF (see Section produces basically the 
same value for the luminosiy density, as all of them are 
in close agreement around the turnover. 

At z = 1, we find that the total rest-frame 
8 /zm luminosity density is (1.4 ± 0.1), (1.3 ± 0.1) and 
(1.4 ± 0.1) X lO^LgMpc^^ for the cases of a double- 
exponential law with cr = 0.36, 0.20 and a Schechter 
function, respectively. This is ^ 4.0 — 4.3 times the cor- 
responding luminosity density at z ~ 0. 

5. THE REST-FRAME 8 Mm LF AT REDSHIFT z - 2 

5.1. The rest-frame 8 fim LF for star-forming galaxies 
at redshift z ~ 2 

The selection of 24 fim galaxies at redshifts z ~ 2 is 
particularly suitable to compute the rest-frame 8/imLF. 
The IR SED of star-forming galaxies is characterised by 
the presence of PAH emission lines from rest-frame wave- 
lengths A = 3.3 through 17/im (Desert, Boulanger & 
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Puget 1990). In particular, one of the main features 
in the SED is the PAH bump around 7.7 and 8.6 //m, 
responsible for a positive selection effect on galaxies at 
z - 1.9 at 24 /im (Caputi et al. 2006a). The light as- 
sociated with this bump produces a substantial fraction 
of the observed 24 /xm output at z 2 (the remaining 
part mainly being produced by AGN). The study of the 
rest- frame 8 ^m LP for star- forming galaxies gives direct 
information on the luminosity distribution of PAH emis- 
sion in IR galaxies. In particular at z ^ 2, it should allow 
us to understand this PAH emission distribution when 
the Universe was only one fourth of its present age. 

We compute the rest- frame 8 /im luminosity {yL%^^) 
of each galaxy in a similar way as for those galaxies at 
0.9<2;<1.1. In this case, the width of the redshift 
bin we consider 1.7 < z < 2.3 implies that the observed 
24 /im maps rest-frame wavelengths 7.2 < Aif < 8.9 /im. 
As we have seen in Section [4. 1[ the k-corrections are usu- 
ally computed based on SED models, which have been 
calibrated using local IR galaxy templates and other ob- 
servables. However, we showed that this can be some- 
what controversial, especially in the PAH-dominated re- 
gion, where different models show important discrepan- 
cies. To compute the k-corrections from 7.2 — 8.9 to 
8 /im , we can avoid relying on any IR SED model by 
directly using measured rest-frame IR spectra of star- 
forming galaxies at redshifts z > 1.5, convolved with the 
24 /im -filter transmission function. These spectra have 
been obtained with the Infrared Spectrograph (IRS) for 
Spitzer (Lutz et al. 2005; Yan et al. 2005). These spectra 
correspond to star-forming ULIRG which are on average 
brighter than those studied here. In spite of that, the 
k-corrections derived within the PAH region from these 
galaxies are expected to be applicable to our galaxies. 
For example, the equivalent widths of PAH lines in the 
Yan et al. (2005) star-forming galaxies are comparable to 
those of other lower luminosity ULIRG. In general, PAH 
line equivalent widths appear to be quite independent of 
the bolometric IR luminosities of star-forming galaxies 
(Peeters et al. 2004; Brandl et al. 2006). 

For the wavelength range considered, the k-correction 
factors derived from empirical spectra vary between k — 
1 (at Arf = 8/im) and k = 1.44 ± 0.36 (at Arf = 8.9 /im). 
These k-corrections are in good agreement with those 
predicted by the Lagache et al. (2004) models. The me- 
dian of the differences is ~ 5% in the considered wave- 
length range (7.2 — 8.9 /im). Thus, the use of empirical 
k-corrections for our rest-frame 8 /im LF at 2 ~ 2 is con- 
sistent with the use of model-dependent k-corrections at 
z=l. 

As at redshift z = 1, we compute the rest- frame 
8/imLF at redshift z ^ 2 alternatively using the 1/Umax 
method and the ML analysis. For the star- forming 
galaxy LF at this redshift, we consider the 132 star- 
forming galaxies with 1.7 < z < 2.3 within our sam- 
ple. The rest-frame 8 /xmLF for star-forming galaxies at 
z ^ 2 computed with the two methods is shown in Fig- 
ure [7] (filled circles for the 1/Umax method and solid and 
dashed lines for the ML analysis). The total comoving 
volume probed at these redshifts is 5.7 x 10^ Mpc~'^. 

For the l/Vmax calculation, we computed the error 
bars taking into account Poisson statistics and the er- 
rors on photometric redshifts and k-corrections through 
Monte-Carlo simulations. We constructed 1000 mock 



catalogues, each one containing 132 galaxies, as the orig- 
inal catalogue. The redshift of each source has been 
allowed to randomly vary following a Gaussian distri- 
bution centred at Zcentrc — z — 0.01 and with a dis- 
persion — 0.06(1 -I- z) (cf. Section \2l3^ . where z is 
the redshift of the source in the original catalogue. The 
redshift of those sources with spectroscopic determina- 
tions have been left unchanged. To include the uncer- 
tainties in the k-corrections, we computed the rest-frame 
8 /im luminosity i^L^'^™ of each galaxy in the mock cat- 
alogue allowing the corresponding k-correction to have 
a random value within the range of its error bar. Once 
more, we recompute the LF with the 1/Vmax method for 
each of the mock catalogues, with exactly the same pro- 
cedure as for the original catalogue. The distribution of 
the LF values in each vLf, ^™ bin determine the error bars 
on our original l/T^nax LF. 

The LF computed with the 1/Umax method which is 
shown in Figure [7] exclusively corresponds to the region 
of z^L^^™ completeness (i^L^^™ > 10" L©). This is es- 
sential to assure that our results are not affected by in- 
completeness effects. 

Also at these redshifts, we analyse the field-to-field 
variations computing the rest-frame 8 /xm LF in the 
GOODS/CDFS and GOODS/HDFN separately. The re- 
sults are shown in the right-hand panel of Figure [9l We 
see that, in spite of the sample variance, the two LF are 
still consistent within the error bars. 

We perform the ML analysis for the combined fields 
in the same way as for galaxies at 0.9 < z < 1.1. Once 
more, we assume that the double-exponential form given 
by eq. (P) with fixed slopes a = 1.2 and l/a^ = 1/(0.36)^ 
can be used to describe the rest-frame 8 /im LF for star- 
forming galaxies at z ~ 2. In this case, the number of 
galaxies is not sufficient to allow us to leave the bright- 
end slope as a free parameter (i.e. the ML algorithm 
does not converge to reasonable values). Also, the adop- 
tion of the same a value as at z ~ is useful to directly 
compare the resulting values of L* and $* at different 
redshifts. Applying the STY method with a double- 
exponential with a — 0.36 to our star-forming galaxies 
at z ~ 2, we obtain that the value of the character- 
istic luminosity which maximizes the ML estimator is 
L* = z/L*8^'^ = (8.3ti:^) X 10^° Lq and the derived nor- 
malization factor is = {9-Otl,]) x 10"" Mpc'^ dex"!. 
The resulting curve for the ML LF at z ~ 2 is indicated 
with a solid line in Figure [71 Once more, the LF ob- 
tained with the ML analysis is in good agreement to that 
computed with the 1/Umax method, confirming that the 
double-exponential form in eq. ([1]) also provides a good 
description of the 8 /imLF also at redshift z ~ 2. 

As for the LF at z = 1, a Schechter function also ap- 
pears to be an alternative suitable law to describe the 
rest- frame 8 /im LF for star- forming galaxies at z ~ 2 
with the ML STY method (dashed line in Figure [T]). 

5.2. Testing the faint-end of the LF through stacking 

analysis 

As we mentioned in Section l473l the faint end slope of 
the IR LF is not well constrained even at z ~ 0, and the 
common procedure is to fix this slope to a given value. 
One could, however, put into question whether the fixed 
slope value we use here (a — 1.2) is realistic to describe 
the faint-end of the IR LF at different redshifts. In the 
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analysis of the IR LF at redshifts < z < 1.2, Le Floc'h 
et al. (2005) concluded that the slope at the faint-end 
could not be much steeper than 1.2, as otherwise the 
faint 24 /im number counts would be overproduced. This 
result has been confirmed by Zheng et al. (2006), us- 
ing the stacking analysis at 24 yum of a large sample of 
0.1 < z < 1 galaxies. The stacking analysis technique 
allows to gain about an order of magnitude in the IR-flux 
detection limit (Dole et al. 2006a; Zheng et al. 2006). 

We do a similar stacking analysis using the Kg < 
21.5 (Vega mag) galaxy sample presented in Caputi et 
al. (2006c) for the GOODS/CDFS. We stack at 24 ^mall 
those galaxies (except AGN) with redshifts 1.7 < z < 2.3 
which are below the completeness limit of the j/L^'^™ lu- 
minosities (i.e. i^i^^™ < 10" Lq at z ~ 2). This in- 
cludes, of course, all those Kg < 21.5 mag galaxies at 
1.7 < z< 2.3 in the GOODS/CDFS which are not iden- 
tified in the S'(24 /xm) > 80 /iJy catalogue for the same 
field. We find that the average 24 /im flux of these stacked 
sources is S'(24/im) = (49.3 ± 1.7) /iJy, which implies an 
average rest-frame 8 yitmluminosity logiQ(i^L^^™) « 10.6. 
To incorporate this stacking point in our differential LF 
expressed per dex unit, we need to estimate the flux - 
and thus the luminosity- range covered by the stack- 
ing sample. Also, wc need to apply a correction factor 
which accounts for the fact that the Kg < 21.5 sam- 
ple loses completeness in identifying 24 ^m galaxies below 
the 5'(24 fj,m) = 80 ^Jy limit. We perform both steps us- 
ing the 24 fim number counts obtained by Papovich et 
al. (2004). These number counts are already corrected 
for incompleteness in the 24 /im detections below the flux 
completeness limits of the Papovich et al. samples. From 
the distribution of these number counts with 24 /im flux, 
we obtain that the average 24 /imflux of our S'(24 /im) < 
80 /iJy sample will be well-reproduced if the stacked 
galaxies span the flux range 30 < 5(24 /im) < 80/iJy. 
On the other hand, from the total number counts within 
this flux range and ignoring the effects of sample vari- 
ance, we can obtain the average identification complete- 
ness produced by the Kg < 21.5 sample. We estimate 
that the Kg < 21.5 sample allows us to identify ^ 79% 
of the 24 /im galaxies with 30 < 5(24 /im) < 80 /iJy. The 
inverse of the completeness fraction gives us the correc- 
tion factor for the LF in the stacking luminosity bin. An 
intrinsic assumption here is that the identification com- 
pleteness derived for this flux range is the same at all 
redshifts, so it can be applied to our 1.7 < z < 2.3 sam- 
ple. This assumption seems to be very plausible (com- 
pare the redshift distributions of the 83 /iJy-limited and 
total 24 /im samples in Figure 3 of Caputi et al. 2006a 
and see also Dole et al. 2006b). 

The resulting stacking point is indicated with a square 
symbol in Figure [71 We note that we only add this point 
to our rest-frame 8 /im LF at z ~ 2 a posteriori, and it 
does not play any role in the ML analysis. The good 
agreement between the stacking analysis point and the 
ML curve confirms that the value fixed for the faint-end 
slope of the 8 /imLF is adequate, and significantly larger 
slopes would not reproduce the average density of faint 
IR galaxies. 

We attempted to do a similar stacking analysis for 
sources at redshifts 0.9 < z < 1.1, in order to test the 
faint-end of the rest- frame 8 /imLF at redshift z = 1. 



However, the stacking at 2A fimoi Kg < 21.5 {Ks < 20.5 
mag) galaxies which are below the 8 /im luminosity com- 
pleteness limit at those redshifts produces an average 
source with flux 5(24 /im) = 16.6 (25.4 /iJy). Unfortu- 
nately, no information on 24 /im number counts are avail- 
able for or below such faint fluxes. This fact prevented 
us to obtain an extension of the rest-frame 8 /im LF at 
z = 1 for faint luminosities. 

5.3. The evolution of the rest-frame 8 /im LF for 
star-forming galaxies from z to z 2 

We can now study the evolution of the rest-frame 
8 /imLF from redshifts z and z = 1 to z ~ 2. Figure 
[8] shows the three LF in a same plot. Different line styles 
in this plot correspond to a double-exponential form with 
(7 = 0.36. As in Section we can characterise the evo- 
lution of L* and $* with redshift. If we consider 



L*(z2 - 2) = L*(zi - 0) X 
$*(z2 - 2) = $*(zi - 0) X 



7L 



1 +Z2 
1 +Zi_ 

1 + Z2^^^ 

1 + Zl 



(6) 



the derived 7l and 75 at strictly zi = 0.2 and Z2 = 1.93 
are 



7L= 3.0 ±0.4 
75 = -2.1 ±0.4. 



(7) 



The obtained 7l value indicates that the strong posi- 
tive luminosity evolution of the rest-frame 8 /im LF con- 
tinues up to redshift z ~ 2. In contrast, the density 
evolution has quite a different trend between z ~ and 
z = 1, and z = 1 and z ~ 2. We showed in Section 
mithat the density of galaxies with ^^L^A^m > iqW.Bj^^ 
dramatically increases from z ~ and z = 1. Between 
and z = 1 and z 2, however, there appears to be a 
significant negative density evolution. If we write 



L*(z2~2) = L*(zi = l) 



$*(Z2'-2) = $*(Z1 = 1)X 



I + Z2 

1 + Zl 

I + Z2 



(8) 



.1 + Zl 

with strictly zi — 1 and Z2 = 1.93, we obtain 
7l= 2.2 ±0.5 

75 = -3.9 ±1.0. (9) 

A negative-density evolution at high (z > 1) redshifts 
have also been found with some of the fittings made for 
the 12 /im LF by Perez-Gonzalez et al. (2005). How- 
ever, these authors conclude that the result of a nega- 
tive density evolution should be taken with caution, as 
it could be produced by incompleteness in the faintest 
luminosity bins. To test this, we repeat the ML analy- 
sis of our rest-frame 8 /im LF by considering only those 
galaxies with 5,y(24/im) > 120/iJy (which is roughly 
equivalent to excluding the faintest luminosity bin in the 
1/Vmax method). In this case, the resulting normaliza- 
tion parameter value implies 75 = — 1.6 ± 0.6 and 
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75 — — 2.7±1.3 for the evolution between z ~ and 
z ~ 2 and between z = 1 and z ~ 2, respectively. We 
conclude, then, that the negative density evolution re- 
sult is not an effect of a plausible incompleteness at the 
faintest luminosities. 

It should be noted that all this analysis is based on 
the validity of the same law to describe the LF at dif- 
ferent redshifts and the values obtained for 7l and 75 
depend on the adopted functional form. A more di- 
rect understanding of the evolution of the rest-frame 
8 fim LF can be achieved by comparing the integrated 
comoving number densities of galaxies above a given lu- 
minosity cut at different redshifts, as those we present 
in Table [5l If we restrict to the most luminous galaxies 
{logiQ{h'Ll'^'^) > 11.5), we find that the number den- 
sity increases with redshift up to z ~ 2. For galaxies 
with logio(j^-^^^™) > remarkably, the number den- 
sity appears to be basically the same at redshifts z = 1 
and z ~ 2. Finally, if we consider those galaxies with 
log]^g(;/L^^™) > 10.5, we observe a clear change of trend 
between z ^ and z = 1, and z = 1 and z 2. While 
the number density of these galaxies increases by a factor 
> 20 between z '--^ and z = 1, then the number den- 
sity at z = 1 decays to half its value by redshift z 2. 
We note that this decrement in intermediate-luminosity 
galaxies is not an effect of the faint-end slope a = 1.2 
we assume for our LF. Inspection of Figure [8] shows that 
only a much higher a value (which would be inconsistent 
with the results of stacking analysis) could make equal 
the number densities of logjQ(z^L® '"") > 10.5 galaxies at 
z = 1 and z ~ 2. 

The rest- frame 8 /xm luminosity density we derive at 
redshift z ~ 2 is 7.5 ± 0.5 (7.6 ± 0.5) x lO^LgMpc-^, ob- 
tained by integrating the double-exponential (Schechter) 
function from the ML analysis. This represents more 
than twice the 8 //m luminosity density at z ~ 0, but 
only half the corresponding luminosity density at z — 1. 

5.4. The total rest-frame 8 /im LF at redshift z ~ 2 

The rest- frame 8 /im LF at z ~ 2 we presented in Sec- 
tion 15.11 has been calculated only taking into account 
the star-forming galaxies in our 24 /itm -selected sample 
at 1.7 < z < 2.3. In this Section, we recompute the 
rest-frame 8/imLF at z ^ 2 for the GOODS fields con- 
sidering all the 161 24 /im -selected star-forming galaxies 
and AGN with 1.7 < z < 2.3. 

We compute the rest-frame 8 /tm luminosities as ex- 
plained in Section [42l To determine the k-corrections for 
the AGN in our sample, we assume that their SED fol- 
low a power-law i/^seu (-^^ith asED < 0). For each 
AGN, we determine the value of asED using its IRAC 
8/imand MIPS 24 /im fluxes. 

The results of the total 8 /tm LF calculated with the 
l/Knax method are indicated with filled diamonds in 
Figure fTOl The error bars take into account Poisson er- 
rors and Monte Carlo simulations on the redshift and 
luminosity catalogues, as explained in Section [5.1l Com- 
paring this total 8 /tm LF with that obtained only for 
star-forming galaxies (Figure ^ , we can see that AGN 
mainly dominate the very bright end. This excess of 
very bright sources suggests that neither the double- 
exponential form given in eq. ([T]) or a Schechter func- 
tion are optimal to describe the bright end of the to- 
tal 8/imLF. At fainter magnitudes, however, the star- 



forming-galaxy and total LF show no significant differ- 
ence, so we can safely assume the same behaviour at the 
faint-end. 

Thus, to compute the total rest-frame 8/imLF with 
the STY method, we consider a combination of an expo- 
nential and a power-law, as follows: 

r $* 1 (J^)'~"x 

I const \ L / 

$(£) = J ^ [-2^^ (1 + i^)] , if i < L* 
[ i^y^"^ , if L > i* 

(10) 

where /3 is the slope at the bright end and the constant 
const = exp{[— l/(2a^)] logiQ(2)} guarantees continuity 
at L = L* . The stacking analysis point (square in Fig- 
ure [ini) is only added a posteriori to check the consis- 
tency of the results. In contrast to the 8 /tmLF for star- 
forming galaxies, the bright end of the total 8/imLF is 
sufficiently well constrained as to allow us to leave f3 as 
a free parameter. At the faint-end, we fix a = 1.2 and 
(7 — 0.36, as in Section [STT] The free-parameter values 
which yield the maximum likelihoood are: /3 = 3.71!q'3, 
L* EE = (2.291°:^^) X 10" Lo and the derived 

normalization is = {3.52to\[l) x lO"'' Mpc'^ dex"!. 
The resulting ML function is plotted with a dashed line 
in Figure [TOl We observe that, while this ML LF is in 
very good agreement with that obtained from the 1/Vmax 
method, the stacking analysis point indicates that the 
faint-end is being under-produced. 

At luminosities 11.0 < logio(i^-Z-^'"°) < 11-4, the 
1/Knax 8/imLF for star-forming and all galaxies are ba- 
sically coincident. However, the slope value a which was 
suitable to describe the former does not seem sufficient 
to explain the faint-end of the total LF. The explana- 
tion for this apparent contradiction is that the values of 
the different free parameters are coupled, and actually 
the definition of faint/bright ends depends on the value 
of L*. In the case of the total rest- frame 8 /tmLF at 
z '--^ 2, the value of the characteristic luminosity L* is 
significantly higher than the ML value of L* for the star- 
forming galaxy LF. We recompute then the STY ML 
estimator for the total LF fixing the slope to a higher 
value a = 1.4. The free-parameter values which yield 
the maximum likelihoood in this case are: (3 = 3.6lto'3, 
L* = uLl^t''^ = (2.34+^:^9) X 10" Lq and the derived 

normalization is = (3.17lo 28) x 10"* Mpc'^ dex"!. 
The ML values of j3 and L* are in agreement with those 
corresponding to a = 1.2, within the error bars. This in- 
dicates the robustness of the determination of the bright 
end and the turnover of the total 8 /tm LF with our sur- 
vey. The resulting ML curve for the case with a = 1.4 
is plotted with a solid line in Figure [TOl This new curve 
appears to be consistent with the stacking analysis point. 

By integrating the obtained STY LF, we can compute 
the 8 /im luminosity density associated with the total IR 
galaxy population at 1.7 < z < 2.3. This luminosity 
density is (9.0 ±0.6) x lO^LgMpc'^, i.e. 2.7 times 
the total 8 /im luminosity density at z 0. Comparing 
the total 8 /im luminosity density at z ^ 2 to that for 
only star-forming galaxies at the same redshift (7.5 ± 
0.5) X 10*" LoMpc^'^, we conclude that AGN have a minor 
contribution to this luminosity density even at high z 
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(- 17% at z - 2). 

6. THE BOLOMETRIC IR LF AT REDSHIFTS z = l AND 

z ~ 2 

6.1. The conversion from vL^^^^^to bolometric L^^j 

6.1.1. A new empirical calibration based on Spitzer 
galaxies 

In Section [5l we studied the rest-frame 8 /imLF at red- 
shift z 2 and its evohition from z 0. The aim of 
this section is to extend this study to the bolometric IR 
(i.e. 5 < A < 1000 fim) LF. The bolometric IR luminos- 
ity of a galaxy is produced by the thermal emission of 
its gas content. In star- forming galaxies, the UV/optical 
radiation produced by young stars heats the interstellar 
dust and the re-processed light is emitted in the IR. For 
this reason, in star- forming galaxies, the bolometric IR 
luminosity allows to obtain a direct and quite unbiased 
estimate of the current star-formation activity. 

Different methods to convert v L^, luminosities into 
bolometric IR luminosities are common in the lit- 
erature. Most of them rely on calibrations made using 
nearby galaxies in IRAS or ISO catalogues (e.g. Chary 
& Elbaz 2001; Elbaz et al. 2002; Takeuchi et al. 2005) 
or on the use of semi-empirical SEDs (e.g. Dale & Helou 
2002; Lagache et al. 2003, 2004; Dale et al. 2005). 
The extrapolation of these ly L^^-L^^ relations to high- 
redshift galaxies can be justified with different recent 
results. For example, Egami et al. (2004) showed that 
composite SED of high-z IR galaxies are well-described 
by local templates. Also, IR-galaxy models which as- 
sume such similarity in the SEDs can fit the 24, 70 
and 160 /xm number counts simultaneously (Lagache et 
al. 2004). Nevertheless, there is always some controversy 
on the large error bars which can be involved in the v L^- 
^hoi conversions applied to high-redshifts. For example. 
Dale et al. (2005) claim that the use of 24 /imdata (rest- 
frame 8/imat z 2) can produce an uncertainty of up 
to a factor five in the derived bolometric IR luminosity 
of z ~ 2 galaxies. 

To explore this issue, Bavouzet et al. (2006) stud- 
ied different v L^-L^^ relations using Spitzer low-to- 
intermediate redshift galaxies. Their sample consists of 
24 /im -selected galaxies with i? < 20 (Vega mag) in three 
different fields, namely the Bootes and the Spitzer First 
Look Survey fields, and the extended CDFS. The se- 
lection criterion of this sample is the detection of each 
galaxy in the IRAC 8 /J,m channel and in all the three 
MIPS bands (i.e. at 24, 70 and 160 ^m). All these galax- 
ies have either spectroscopic or C0MB017 photometric 
redshifts and span the redshift range 0.0 < z < 0.6. AGN 
have been removed from their sample. 

To measure the bolometric IR luminosity L^^j of each 
galaxy at redshift z, Bavouzet et al. (2006) used the 8 
through 160 /im fluxes. To integrate the corresponding 
empirical SED in each case, they summed up the areas 
below contiguous rectangles centred at rest-frame wave- 
lengths 8/(l-Hz), 24/(l-hz), 70/(l-Fz) and 160 A<m/(l-hz). 
At longer wavelengths, they approximated the SED be- 
yond A > [160 H- (160 - 70)/2]/(l + z) = 205 /im/(l + z) 
with a triangle of slope —4. This slope is consistent 
with the modified black-body emission in the far-IR pro- 
duced by big grains of dust thermalised at a temperature 
T ~ 15 - 20K (Draine & Lee 1984; Contursi et al. 2001). 



In fact, Bavouzet et al. (2006) found that the use of any 
slope between —3.5 and —4.5 produced variations < 1% 
on the computed bolometric luminosities. It is impor- 
tant to emphasize that the measurements of bolometric 
IR luminosities made by Bavouzet et al. are purely based 
on Spitzer data and are completely independent of any 
model template. 

The resulting versus rest-frame z/L^ '""relation for 
the Bavouzet et al. sample is shown in Figure [Tl] (cross- 
like symbols). In this work, however, we restrict the 
analysis only to those 93 galaxies in the Bavouzet et 
al. sample which have z/L^^™> 10^° L© and signal-to- 
noise S/N > 3 ratio in all the MIPS bands. The rest- 
frame 8 /im luminosities have been obtained applying k- 
corrections which do depend on different SED models 
(Chary & Elbaz 2001; Elbaz et al. 2002; Lagache et 
al. 2004). The z/L^Mm_LiR relation for these galaxies 
can be fitted with the following law (dashed line in Fig- 
ure E]): 

L™i. = 1.91 X [^.L,(8/im)]l■0^ (11) 

with i^Lij (8 nm) and expressed in units of Lq . The 

la dispersion for this relation is 55%. This formula is 
directly applicable in all the redshift range 0.0 < z < 0.6. 
We refer the reader to the Bavouzet et al. (2006) pa- 
per for a generalised version of this formula including 
vLl^''^< lOi°L0 galaxies. 

To assess whether the formula displayed in eq. 
(fTT|) could also be suitable to perform the i^L®^™- 
L^^^ conversion for higher redshift galaxies, Bavouzet et 
al. (2006) used the 24 /im -selected galaxy samples in the 
GOODS/CDFS and HDFN (the same samples we use 
in this work). Of course, the bolometric luminosity of 
the vast majority of z > 1 galaxies cannot be empiri- 
cally measured, as they are below the confusion limits of 
the Spitzer /MIVS images at 70 and 160 /im. However, 
the average far-IR flux produced by these high-redshift 
sources can be recovered through stacking analysis (Dole 
et al. 2006a). 

Bavouzet et al. (2006) stacked all those 24 /im sources 
in the GOODS fields that he at redshifts 1.3 < z < 2.3, 
with a median redshift z « 1.68. The resulting (i/L^'^'"; 
^boi ) value obtained with the stacking analysis is indi- 
cated with a filled circle in Figure [11] The value 
for this point is corrected for the subestimation of the 
far-IR flux which is produced on high redshift sources 
by using the triangle-approximation method explained 
above. This correction is about 10 — 15%. The lo- 
cus occupied by the high-reshift stacked sources in the 
^j^^^im jiK^ diagram strongly suggests that the relation 
described by eq. (fTTjl is also valid to link the S/^mand 
bolometric IR luminosities of IR galaxies at 1.3 < z < 
2.3. 

Thus, in this work we make use of the Bavouzet et 
al. relation described by eq. pT|) to convert the rest- 
frame 8 /im of our star- forming galaxies into bolometric 
IR luminosities. We use these resulting bolometric IR 
luminosities to compute the corresponding LF for star- 
forming galaxies at z = 1 and z ~ 2 that we present in 
Section [621 As we explain in that section, the 55% dis- 
persion found for the relation largely domi- 
nates the error budget of the bolometric IR LF. 

As a final comment, we would like to discuss why 
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the relation obtained by Bavouzet et al. (2006) pre- 
dicts a significantly smaller dispersion on the values of 
bolometric IR luminosities obtained from rest-frame 
8 /xm fluxes than that predicted by Dale et al. (2005). 
Firstly, the Dale et al. (2005) sample includes nearby 
galaxies of very different nature, and they even make 
separate studies of different IR regions within a same IR 
galaxy. Thus, because of its selection, it is expected that 
the Dale et al. sample displays a larger variety of IR 
properties than the Bavouzet et al. (2006) sample. Fur- 
thermore, to extrapolate their conclusions to high red- 
shifts. Dale et al. (2005) make use of the complete set 
of Dale & Helou (2002) models. However, the majority 
( ^ 75%) of their wide range of observed SEDs only cor- 
respond to roughly half of these models (cf. figures in 
Dale et al. 2005). The Bavouzet et al. sample has been 
selected with a more homogeneous criterion and includes 
galaxies up to intermediate {z ss 0.6) redshifts. Thus, 
these galaxies are more likely representative of the typi- 
cal galaxies selected in IR surveys. A thorough discussion 
of this issue is presented in the Bavouzet et al. (2006) pa- 
per. 

6.1.2. Comparison between different i/L^''™- 
LJjqj calibrations 

Several different laws to convert i/L^'^™into bolometric 
IR luminosities are of common use in the literature. 
We analyse here how these different calibrations compare 
to the relation empirically derived from Spitzer galaxies 
by Bavouzet et al. (2006). 

Figure [12] shows the bolometric IR versus 
i^L^ ^""luminosity relations (left-hand panel) and the de- 
rived conversion factors /(^^^^'"(right-hand panel), 
both versus I'L^'^™, as obtained using different calibra- 
tions or derived from different IR SED models. The 
thick solid line shows the empirical relation obtained 
by Bavouzet et al. (2006). The thick dashed and dot- 
ted lines correspond to the relations derived using the 
Lagache et al. (2004) and the Chary & Elbaz (2001) 
- Elbaz et al. (2002) templates, respectively. To ob- 
tain these relations, we convolve all these templates 
with the transmission function of the MIPS 24 /im filter. 
We find that the Lagache et al. model predicts a 
i,]^Sfj.m_ j^m^ relation which is in close agreement with 
the Bavouzet et al. (2006) empirical calibration over 
all 8 /im luminosities. The Chary & Elbaz templates, 
on the contrary, appear to over-produce the i/L^'^™- 
^hoi conversion. The differences with the Bavouzet et 
al. calibration are only within a factor ~ 2 for galaxies 
with z/L^''™ < 10^^ L0, but become dramatically larger 
at higher luminosities. 

Previous comparisons of the v L^^-L^^ relations pre- 
dicted by different models have not detected such dra- 
matic differences (see e.g. Le Floc'h et al. 2005). These 
previous comparisons analysed longer rest-frame wave- 
lengths, beyond the PAH-dominated region in the SEDs. 
The comparison we present here is made in the most 
critic SED region, where different models show the 
largest discrepancies (cf. also Figure [5|). From this com- 
parison, we find that the use of the Chary & Elbaz tem- 
plates to convert i^L^'^^into luminosities leads to 
significantly over-produced bolometric IR luminosity val- 
ues for galaxies with i/Lf,'^™ > 10^^ Lq. 



In Figure fT2l we also show the i^L*^™to derived 
from the Wu et al. (2005) formulae (thin dashed lines), 
which link 8 /im luminosities and star formation rates. 
The bolometric IR luminosities have been recovered us- 
ing SFR = 1.72 X 10-^° Lin (Kennicutt 1998). Fi- 
nally, the thin dotted-dashed line shows the relation used 
in Reddy et al. (2006a). In this latter relation, the 
8 /im luminosities refer to the convolution in the wave- 
length range ~ 5 — 8.5 /xm, which is somewhat different 
from the transmission windows of the MIPS 24 /im filter 
(- 6.6 - 9.4 /im at z ~ 2) or the IRAC 8/imfilter 
6.5 — 9.5 /xm; Fazio et al. 2004). Once corrected for 
this difference, the Reddy et al. (2006a) relation becomes 
closer to the Bavouzet et al. (2006) Spitzer calibration. 

In this work use the new Spitzer-hased calibration 
given by eq. (fTT|) to convert i/L^ ^""luminosities into bolo- 
metric IR luminosities . After computing the bolo- 
metric IR LF, we analyse the contribution of LIRG and 
ULIRG to the total number and luminosity densities of 
IR galaxies at different redshifts. We warn the reader, 
however, on the implications of the differences between 
the vLf,'^™- conversions shown in Figure fT2l For ex- 
ample, the Chary & Elbaz conversion classifies as ULIRG 
to those sources with i/L^^™ ^ 8 x 10^° L©, while the 
Bavouzet et al. relation implies that only galaxies with 
^2,8 Min > _ 1.2 X 10" Lq are ULIRG. These differ- 
ences should be kept in mind when comparing differ- 
ent results from the literature, where different conversion 
laws are used. 

6.2. The bolometric IR LF for star-forming galaxies 
and its evolution to redshift z ^ 2 

As we have seen in Section 16.1.21 some calibrations 
used in the literature to convert 8 /im into bolometric IR 
luminosities are quite discrepant with the empirical cal- 
ibration obtained from Spitzer galaxies. Thus, to prop- 
erly compare the bolometric IR LF at different redshifts, 
we convert the different 8 fim LF using the Bavouzet et 
al. (2006) relation shown in eq. (jlip . The results are 
shown in Figure [T51 

Firstly, we transform the Huang et al. (2006) 8 /imLF 
at z ~ and compute the corresponding minimum x^- 
fitting, using the functional form given in eq. For 
the bolometric IR LF at z ~ 0, we obtain the following 
best-fit parameter values: a = 0.39 ± 0.01, Ljp — (4.0 ± 
0.2) X 10i°Lo and $* = (5.4 ±0.1) x lO^^ ^pc-s dex"!. 
The resulting best-fit curve to the bolometric IR LF at 
z ^ is shown with a dotted line in Figure [131 

The best-fit value we find for the slope at the bright- 
end at z 0, i.e. a = 0.39, is very similar to the value ob- 
tained for the bright-end slope of the rest-frame 8 /im LF 
(cr — 0.36) at the same redshift. This similarity is due to 
the fact that the conversion is quasi-linear. 

At redshifts z = 1 and z 2, we compute the bolo- 
metric IR luminosities L^^ of all our star-forming galax- 
ies in the relevant redshift ranges by transforming their 
rest- frame 8 /xm luminosities z^L^^'"using eq. ([T]). We 
then obtain the bolometric IR LF using, alternatively, 
the l/Vinax method and the ML STY analysis. 

The upward-pointing triangles and circles in Figure 
[T3l show the bolometric IR LF at z = 1 and z ^ 2, 
respectively, both computed with the 1/Vmax method. 
These LF are only shown in the bins of completeness 
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in Lj^^j luminosities, given the flux limits of our sample 
and the redshifts corresponding to each case. As for the 
rest- frame 8/j,mLF, we applied small correction factors 
to account for the 5 — 6% identification incompleteness of 
the S'(24/im) > 80 /iJy galaxy sample. For both LF, the 
error bars have been determined through Monte Carlo 
simulations, in a similar way as described in Section f5.ll 
The mock catalogues generated in the simulations take 
into account the error bars in the redshift determina- 
tions, in the case of photometric redshifts. However, in 
the case of the bolometric luminosities, the error budget 
is mainly dominated by the uncertainty associated with 
the vL^^^'°^-L^Y conversion. To take into account this 
error, we assign to each galaxy in the mock catalogues a 
random bolometric IR luminosity. This random luminos- 
ity L{j^j is taken from a Gaussian distribution centred at 
the value given by eq. pT|) for the corresponding galaxy 
and with a 55% dispersion. The re-computation of the 
LF with the l/Knax method on all the mock catalogues 
allows us to determine the error bars on the original LF 
calculation. 

The dotted-dashed and solid lines in Figure [T3l indicate 
the bolometric IR LF at z = 1 and z ~ 2, respectively, 
obtained with the ML analysis. We computed the bolo- 
metric IR LF using the STY method, assuming the func- 
tional form described in eq. Ilj . The faint and bright-end 
slope values have been fixed to the z ~ values, i.e. a — 
1.2 and a = 0.39, respectively. At z = 1, we obtain that 
the value of the characteristic luminosity which yields the 
ML is = (2.5ti^;3) x 10" ^0- The corresponding nor- 
malization factor is $* = {^-Qtl'l) x lO'^ Mpc'^ dex"!. 
At z 2, the ML characteristic luminosity is L^^ = 
(6.3lo g) X 10^^ L0 and the corresponding normaliza- 
tion factor is $* = (9.2t^;^) x IQ-* Mpc'^ dex"! (cf. 
Table [7|). The error bars on ijp include the uncer- 
tainty produced by the 55% dispersion in the 
^boi relation, incorporated through the mock catalogues 
described above. Consistently with the results obtained 
in Sections 14. 31 and 1 5.11 the LF independently calculated 
with the l/Knax method and the ML STY technique are 
in very good agreement. 

Using also the vLl'^^-L^i relation given in eq. 
(jlip . we compute the corresponding contribution of 
the stacked galaxies at z ~ 2 which are below the 
completeness limit of the sample, to the bolometric 
IR LF. Once more, the stacking analysis point appears 
in very good agreement with the extrapolation given by 
the ML analysis at the faint-end of the LF. 

Given the quasi- linearity of the vLf,'^'^- 
^hoi conversion, the evolution we find for the bolometric 
IR LF from z ~ to z ~ 2 is very similar to the 
evolution observed for the rest-frame 8/xmLF. For the 
bolometric IR LF, this implies: 

• The number density of galaxies with 
-^boi ^ 10^^ L0 substantially increases from 
the local Universe to z = 1 (cf. Table!]). This 
confirms the increasing importance of the LIRG 
and ULIRG populations between these redshifts 
(cf. e.g. Le Floc'h et al. 2005). 

• Surprisingly, at z ~ 2, the number density of star- 
forming ULIRG (i.e. sources with > 10^^ Lq) 



is only slightly larger than at z = 1. This result 
is the combination of several factors: firstly, the 
exclusion of AGN in this analysis produces a rela- 
tively low density of ULIRG at z ~ 2, as we have 
seen in Section 15.41 that AGN dominate the bright 
end of the IR LF. Secondly, the use of the i^L^^'"- 
to-iJ,^j conversion given in eq. (fTT|) . which, m com- 
parison to the Chary & Elbaz templates that are 
of common use in the literature, produces ULIRG 
only from larger zyL^^^luminosities (cf. Figurc [T2)) . 

• The number density of LIRG (i.e. sources with 
10^^ < < 10^^ L0) appears to be smaller at 

z 2 than at z = 1. Although the limits of our 
survey do not allow us to directly observe LIRG at 
z ~ 2, the ML analysis suggests this results, which 
is in turn validated through the stacking analysis 
of z ~ 2 Ks-hand galaxies. 

Thus, the ratio between the number densities of star- 
forming ULIRG and LIRG increases from z = 1 to z ~ 2. 
However, within our sample and given our star-forming 
galaxy/ AGN separation, this effect appears to be mainly 
produced by a decrement in the density of LIRG by 
z ~ 2, rather than a significant increment in the density 
of star- forming ULIRG. If our AGN separation criterion 
were excluding galaxies whose bolometric IR emission 
is actually dominated by star-formation, then the rela- 
tive importance in the number density of star-forming 
ULIRG would be, of course, even larger at z ~ 2. 

We note that the decrement we find in the number den- 
sity of LIRG between z = 1 and z ~ 2 is not influenced 
at all by the AGN separation criterion. 

6.3. Comparison with other works 

As we have seen in Section 16.1.21 many different 
recipes are used in the literature to convert v into 
-^boi luminosities. And even different conversions made 
from a same wavelength (in particular, rest-frame 8 /im ) 
may lead to non-negligible discrepancies in the derived 
luminosities. In spite of these differences, it is still 
instructive to compare the results of different bolometric 
IR LF calculations. 

Figure [14] compares the bolometric IR LF obtained in 
this work with those derived by other authors, at dif- 
ferent redshifts. In the left-hand panel, we show the lo- 
cal bolometric IR LF computed from the IRAS revised 
galaxy sample (Sanders et al. 2003; diamond-like sym- 
bols) and the bolometric IR LF derived in this work from 
the Huang et al. (2006) rest-frame 8/j,mLF a.t z ^ 0.2. 
The difference between the two is mainly due to a real 
evolution between z = and z ~ 0.2. In the same panel, 
we also compare our bolometric IR LF at z = 1 with 
that obtained by Le Floc'h et al. (2005) at z = 0.9. We 
observe that both LF are in good agreement, taking into 
account the error bars and the evolution expected be- 
tween these redshifts (cf. Le Floc'h et al. 2005). 

In the right-hand panel of Figure [14] we show our bolo- 
metric IR LF at redshift z ~ 2, compared to that derived 
from Perez- Gonzalez et al. (2005) at a similar redshift 
and that computed from radio-detected sub-millimetre 
galaxies at z ~ 2.5 (Chapman et al. 2005). 

The bolometric IR LF derived from Perez- Gonzalez 
et al. (2005; asterisks in Figure [14)) has been obtained 
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by converting their rest-frame 12 fim LF at z ^ 2, us- 
ing the same recipe adopted by these authors to ob- 
tain bolometric IR luminosity densities (cf. eq.(l) in 
their paper). This conversion corresponds to the Chary 
& Elbaz (2001) vLl^ ^''^-L^^^ ^orIml\a,. Our bolometric 
IR LF at z ~ 2 is in agreement, within the error bars, 
with that derived from Perez- Gonzalez et al. (2005) at 
luminosities Lj^^j < 10^^'^ L©. At brighter luminosi- 
ties, however, the two LF present significant discrepan- 
cies. The differences between the two are produced by 
two factors: 1) the AGN exclusion: Perez-Gonzalez et 
al. (2005) only exclude the most extreme cases of AGN, 
while here we adopt a more extensive separation crite- 
rion; 2) the different v L^^-L^^ conversions: as we have 
seen in Section [5.1.2[ the most drastic differences between 
the empirical Spitzer-hased conversion we use in this 
work and that derived from the Chary & Elbaz (2001) 
and Elbaz et al. (2002) templates occur at luminosities 
-^boi ~ 10^^ L0. This comparison illustrates the impact 
of using different v L^-L^^ relations, especially at high 
redshifts, where the most luminous IR galaxies are dom- 
inant. 

The bolometric IR luminosities derived from radio- 
detected sub-millimetre galaxies only can trace the very 
bright-end of the bolometric IR LF. The diamond- like 
symbols in the right panel of Figure [14] correspond to 
the sub-millimetre-derived bolometric IR LF at z 2.5, 
as obtained by Chapman et al. (2005). This LF does 
not exclude AGN and quickly loses completeness at 
-^boi. ~ 10^^ L0. Taking into account these facts and 
the differences in redshift, we find that the Chapman et 
al. (2005) bolometric IR LF at z ^ 2.5 is consistent with 
our LF based on 24 /im -selected galaxies at z ~ 2. 

6.4. The evolution of the bolometric IR luminosity 
density 

One of the final aims of computing the bolometric IR 
LF is to obtain an estimate of the IR luminosity den- 
sity (in our case associated with star-forming galaxies) 
at a given look-back time. Previous works agree in a 
strong evolution of the IR luminosity density from the lo- 
cal Universe up to redshift z ~ 1 (e.g. Flores et al. 1999; 
Gispert, Lagache & Puget 2000; Pozzi et al. 2004; Le 
Floc'h et al. 2005). At higher redshifts, the situation is 
less clear as only recently IR facilities are allowing to put 
constraints on the IR Universe at z > 1. 

Given the discrepancies existing between different 
recipes to obtain bolometric IR luminosities (cf. Section 
I6.1.2p . we need to use the bolometric IR LF obtained 
with the same conversion at different redshifts, in order 
to properly compute the evolution of the IR luminosity 
density. 

Figure [15] shows the evolution of the comoving IR 
luminosity density as a function of redshift. Our de- 
terminations of the IR luminosity density at z = 1 
and z ^ 2 (strictly z = 1.93) are indicated with a 
filled upward-pointing triangle and circle, respectively: 
nin{z = 1) = (1.2 ± 0.2) X lO^LgMpc-^ and n^iz ~ 
2) = (6.6lit;?) X lO^LgMpc-^. We obtain the values of 
these luminosity densites by integrating our respective 
bolometric IR LF obtained with the ML likelihood anal- 
ysis, weighted with the luminosity values. The error bars 
are determined by the extreme cases of LF produced by 



the error bars on L^^. 

The cross- like symbol in Figure [15] represents the bolo- 
metric IR luminosity density at z ~ (strictly z = 0.2), 
as obtained from the bolometric IR LF derived from the 
Huang et al. (2006) 8^mLF: nm{z - 0) = (2.5 ± 0.2) x 
10*LoMpc"3. 

The thick solid line in Figure [T^] interpolates the evo- 
lution of the total bolometric IR luminosity density be- 
tween redshifts zi ~ to Z2 = 1 and zi = 1 to Z2 ~ 2, 
assuming this evolution follows a [(1 -t- Z2)/(l -I- zi)]^ 
law. Between redshifts z ~ and z = 1, we find that 
the total bolometric IR luminosity density increases as 
[(1 + Z2)/(l + zi)]3-i±0-3 (where zi = 0.2 and Z2 = 1.0). 
This evolution is somewhat slower than that obtained by 
Le Floc'h et al. (2005), who found [(1 + Z2)/(l + zi)f-^ 
between zi = and Z2 = 1. The bolometric IR luminos- 
ity density at z = 1 determined by Le Floc'h et al. (2005; 
right-hand pointing triangle in Figure [T5)) is actually very 
close to the value we determine here. The difference ap- 
pears to be mainly produced in the IR luminosity density 
at low redshifts: there has been a significant evolution of 
the IR LF between redshifts z = and z ~ 0.2. 

Other symbols in Figure[T5]refer to different bolometric 
IR luminosity density estimations derived from different 
datasets: radio (Haarsmaet al. 2000; downward-pointing 
triangle), sub- millimetre (Barger et al. 2000; diamond- 
like symbol) and the different fits made on mid-IR data 
by Perez-Gonzalez et al. (2005; asterisks). Our determi- 
nations of the IR luminosity densities are in good agree- 
ment with most of these previous works within the error 
bars. Our results exclude, however, the highest of the 
three estimations made by Perez-Gonzalez et al. (2005) 
at z > 1. 

Finally, in Figure [15] we show the relative contribu- 
tions of the LIRG and ULIRG populations to the to- 
tal IR luminosity density, as a function of redshift. At 
z ~ 0, (28^20)*^ of bolometric IR luminosity density 
is contained in LIRG and < 1% in ULIRG. At z = 1, 
we find that LIRG and ULIRG contribute (6117)% and 

(16^12)%, respectively, to the total IR luminosity den- 
sity, in agreement with Le Floc'h et al. (2005) within 
the error bars. By z ^ 2, the contribution of LIRG and 
ULIRG become iA7tn)% and (42t^^)% of the total bud- 
get, respectively. 

Using the Kennicutt formula SFR = 1.72 x 10"^° Lir, 
we can convert the bolometric IR luminosity densities 
into star- formation rate densities at different redshifts. 
At z = 1 and z 2, respectively, fiiR = (1.2 ± 
0.2) X 10^ and {G-GtH) x IO^LqMpc-^ translate into 
star-formation rate densities (5sfr ~ (0.20 ± 0.03) and 
(0.11 ±0.02) Mgyr"! Mpc'^ (assuming a Salpeter initial 
mass function over stellar masses M = (0.1 — 100) Mq). 
In Section [T] we make use of our current knowledge on 
stellar mass density evolution to discuss why these de- 
rived star-formation rate densities could not be much 
higher than this value at redshifts 1 < z < 3. 

7. DISCUSSION 

If the IR LF for star-forming galaxies follows a unique 
law from the local Universe to high redshifts, then the 
results of our LF determination will imply that there 
is a negative evolution in the overall number density of 
IR star-forming galaxies between z '--^ and z 2. We 
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showed here the vahdity of a universal law to describe the 
IR LF at intermediate and bright luminosities, at differ- 
ent redshifts. Of course, one could argue that the faint- 
end of this LF is not sufhciently well constrained as to 
determine the number density of low-luminosity objects. 
Although a direct probe of the faint-end of the IR LF will 
require the capabilities of next-generation telescopes as 
the James Webb Space Telescope (JWST), the stacking 
analysis of galaxies below the limits of our 24 /im survey 
appears to support our conclusion. The result of stacking 
analysis suggests that the faint-end slope of the IR LF at 
z ~ 2 cannot be much higher than the value we consid- 
ered here (and those usually considered in the literature 
at different redshifts). 

In fact, an analogous situation is observed at other 
wavelengths. For example, Caputi et al. (2006c) deter- 
mined the evolution of the rest-frame i^s-band LF from 
z = to z ~ 2.5. The depth of their survey {Kg < 21.5 
Vega mag) allowed them to properly constrain this LF 
down to more than a magnitude below the turnover M* 
at z = 2. These authors found that a Schechter func- 
tion with a same fixed slope is suitable to describe the 
Ks-hand LF from the local Universe to high redshifts, 
within the limits of their survey. In this case, the ML 
analysis (which is in good agreement with the LF com- 
puted with the 1/Knax method) also indicates there is 
a negative density evolution of this LF with increasing 
redshift. 

The similarities between the evolutions of the Kg-hand 
and 8 /imLF should not come as a surprise. The bright- 
end of the mid-IR LF at z ~ 2 is mostly populated by 
massive M > 10" Mq galaxies (Caputi et al. 2006a). At 
redshift z ~ 1, the mid-lR LF is dominated by LIRG, 
the majority of which are characterised by intermediate 
~ 10^° - 10" Mq stellar masses (Hammer et al. 2005; 
Caputi et al. 2006a). Thus, the evolution in the number 
density of mid-IR galaxies above a given luminosity cut 
is related to the global evolution of galaxies above a given 
mass cut. 

It should be clear that the aim of this discussion is to 
show how the results wc find in this work are perfectly 
consistent with other observational evidence of galaxy 
evolution. This does not exclude, however, that the ulti- 
mate conclusion on the faint-ends of the Kg and IR LF 
will only be achieved in the light of future extremely deep 
surveys. 

As we mentioned in Section 16.41 the IR luminos- 
ity density associated with star-forming galaxies at 
z 2 implies a star- format ion rate density (0.11 ± 
O.O2)M0yr"iMpc~^ (Kennicutt 1998). Let us assume 
that this has been the average star-formation rate den- 
sity between redshifts z = 1 and z = 3. In our assumed 
cosmology, the elapsed time between these redshifts is ~ 
3.6 Gyr. The stellar mass density formed during this pe- 
riod of time would be nearly (4.0 ±0.7) x 10® Mg Mpc"^. 
If we consider that the fraction of material recycled 
through stellar winds and supernovae could be roughly 
50%, then the resulting mass density locked in stars 
would grow by - (2.0 ± 0.4) X 10* Mq Mpc'^ between 
z = 3 and z = 1. This is actually the growth of the stel- 
lar mass density that has been measured from different 
near-IR surveys at these redshifts (see Caputi et al. 2006c 
and references therein). This result also shows that, un- 
less the recycled fraction of material into the intestellar 



medium is much larger than 50% between redshifts z = 1 
and z = 3, then the average star-formation rate density 
along this period cannot very much exceed the value we 
find in this work Ssfr « (0.11 ± 0.02) Mq yr'^ Mpc"^ at 
z ~ 2. Much higher star- formation rate densities only 
could be explained if a high fraction of the new formed 
stars were very massive, in which case they would not 
basically contribute to the final stellar mass of the host 
galaxies. 

Considering a star formation rate density (5sfr = 
(0.11±0.02) Mq yr~i Mpc"^ strictly in the redshift range 
1.7 < z < 2.3 and assuming again a recycled frac- 
tion of 50%, we derive that the total stellar mass den- 
sity produced in this redshift interval is (1.8 ± 0.3) x 
IO^MqMpc"^. This is nearly 4% of the total stellar 
mass density assembled at z = (i.e. (4.9 ± 0.1) x 
10* Mq Mpc~'^, as obtained by integrating the local stel- 
lar mass function of e.g. Cole et al. 2001). In the redshift 
interval 0.9 < z < 1.1, our measured star formation rate 
density is Ssfr = (0.20 ± 0.03) Mq yr"! Mpc'^. With a 
recycled fraction of 50%, this implies a growth in stellar 
mass density of (8.0 ±1.2) x lO^MgMpc"^. Thus, more 
than 15% of the present-day stellar mass density is being 
created in IR galaxies during the time elapsed between 
redshifts z = 0.9 and z = 1.1 (i.e. ~0.8 Gyr). 

We found in this work that the number densities of 
ULIRG associated with star formation are very similar 
at redshifts z = 1 and ^ 2. This suggests that the 
physical mechanism responsible for galaxies to enter a 
star-forming ULIRG phase is similarly efficient at these 
two redshifts. This result imposes strong constraints on 
IR-galaxy synthesis models. The origin of the ULIRG 
phase is usually associated with advanced gas-rich merg- 
ers (Sanders & Mirabel 1996). Thus, this phenomenon 
had to be comparably common for the production of 
powerful star-forming systems at redshifts z = 1 and 
2. 

8. SUMMARY AND CONCLUSIONS 

In this work we have presented the IR LF of 24 /im - 
selected Spitzer galaxies at redshifts z = 1 and z ^ 2 in 
the GOODS fields. At z ~ 2, we separately studied the 
LF for star-forming galaxies only and the total 8 /im LF 
for star-forming galaxies and AGN. We then used a new 
calibration based on Spitzer star-forming galaxies to con- 
vert the rest- frame 8 /im into bolometric IR luminosities 
of the star-forming galaxies in our sample. This allowed 
us to compute the bolometric IR LF and obtain an esti- 
mate of the IR luminosity densities at z — 1 and z ~ 2. 

We found that the rest-frame 8 /im LF for star- forming 
galaxies at z = 1 and z ~ 2 is well described by a double- 
exponential law which has evolved from z ~ 0. Between 
z ~ and z = 1, there is a strong luminosity evolu- 
tion and the number density of log]^Q(;/L*'"") > 10.5 in- 
creases by a factor > 20. The characteristic luminosity 
L* of the rest-frame 8 /xmLF continues increasing up to 
redshift z ~ 2, but, at this redshift, the number density 
of log]^Q(;/L*'"") > 10.5 galaxies is smaller than the den- 
sity at z = 1. This certainly does not mean that the 
contribution of IR galaxies has been less important at 
high redshifts. The rest-frame 8 /im luminosity density 
at z ^ 2 is still ^ 2.3 times larger than the correspond- 
ing luminosity density at z ~ 0, but only half the value 
at z = 1. 



16 



K. I. Caputi et al. 



At z 2, the inclusion of AGN mainly affects the 
bright end of the IR LF. The bright end of the total 
rest- frame 8 /um LF for star-forming galaxies and AGN is 
correctly reproduced by a power-law which accounts for 
the excess of bright sources. AGN only produce ~ 17% 
of the total rest-frame 8 /lim luminosity density at 2; ~ 2. 

The quasi-linear relation between rest- frame 8 /im and 
bolometric IR luminosities for star-forming galaxies 
makes that the bolometric IR LF is well-described by 
a similar law as the rest-frame 8 /zm LF at the same red- 
shift. The characteristic luminosity L^^ of the bolomet- 
ric IR LF for star-forming galaxies at z 2 is close 
to ~ 10"'^^ Lq, i.e. the limiting luminosity between the 
LIRGs and ULIRGs. As the luminosity density is mainly 
governed by the turnover of the LF, the value of L^^i 
results in roughly similar contributions of LIRGs and 
ULIRGs to the IR luminosity density. These two popu- 
lations altogether account for ~ 90% of the total IR lu- 
minosity density associated with star formation at 2; ~ 2. 

Finally, we discussed the possibility that the total IR 
luminosity and corresponding star-formation rate den- 
sity estimated in this work could have been significantly 
different at any redshift between z — 1 and z = 3. Con- 
straints from near-IR surveys suggest that the stellar 
mass density built up by galaxies at this epoch would be 
in contradiction with average star-formation rate densi- 
ties much larger than our estimated value (unless a much 
higher proportion of very massive stars were created in 
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TABLE 1 
The rest-frame 8 /xm LF for 

STAR-FORMING GALAXIES AT 2 = 1 
OBTAINED WITH THE 1/Vmax METHOD. 

logloCi^Lg^-") logio$(Mpc-3dex-^) 

+0.06 
-0.08 
+0.07 
-0.07 
+0.09 
-0.10 
+0.12 
-0.15 
+0.25 
-0.54 



10.60 


-2.55; 


10.80 


-2.66; 


11.00 


-2.97; 


11.20 


-3.30; 


11.40 


-4.12" 



TABLE 2 

Parameter values characterising the rest-frame 8/.imLF for star-forming galaxies 

AT 2 = 1. 

Functional form a a uLl^ ^"'^^ (Lq) $* (Mpc-^dex"!) 

Double-exp (eq. [TJ 1.2 (fixed) 0.36 (fixed) (3.55^0,11) x 10^° (3.95lo jg) X lO'^ 

Double-exp (eq. [TJ 1.2 (fixed) 0.20to:ln (free) (1.10to;g|) ^ 1°" (2-54lo 3°) X 10"^ 

Schechter (eq. [I 1.2 (fixed) — (7.2^55;^) X 10^1^ (3.88lo 4i) X lO'^ 



TABLE 3 
The rest-frame 8 /^m LF for 

STAR-FORMING GALAXIES AT 2 ~ 2 
OBTAINED WITH THE l/Vmax METHOD. 

Iogio(t/Lg^'") Iogio$(Mpc-3dex-i) 
11.09 -3.34+1-^1 



11.29 



1+0.09 



11.49 -3.oa_o.„ 
11.69 -4.58+"'^^ 



+0.18 
0.13 
0.29 
0.38 



TABLE 4 

Parameter values characterising the rest-frame 8 fim LF for, star-forming 

GALAXIES AT 2 ~ 2. 



Functional form 


a a 


^iD'^-" (Lo) 


(Mpc-^dex-i) 


Double-exp (eq. [TJ 
Schechter (eq. [Sj 


1.2 (fixed) 0.36 (fixed) 
1.2 (fixed) — 


{,8.3+\i) X lOio 
{lS2t°il) X 10" 


(,9.0+f l) X 10-4 
{9.3tl:l) X 10-4 
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TABLE 5 

Number densities of galaxies with rest-frame i/LJ'""above different luminosity cuts 
AT different redshifts. These number densities have been obtained by integrating 

THE FUNCTIONAL FORM APPEARING IN THE SECOND COLUMN AND ARE EXPRESSED IN UNITS OF 
Mpc^^. DE STANDS FOR DOUBLE- EXPONENTIAL. 



Redshift 


Functional form 


log^oi'^Lt'^n > 10 


.5 


> 11.0 




> 11.5 




2 ~ 


DE (cr = 0.36) 


(4.8 ±0.4) X IQ- 


-5 


(6.7 ±0.9) X 10" 


-7 


(1.4 ±0.3) X 10- 


-9 


2=1 


DE (o- = 0.36) 


(1.1 ±0.1) X 10- 


-3 


(1.8 ±0.3) X 10- 


-4 


(6.7 ±2.0) X 10- 


-6 




DE (o- = 0.20) 


(1.1 ±0.1) X IQ- 




(1.7 ±0.3) X 10- 


-4 


(2.1 ± 1.0) X 10- 


-6 




Schechter 


(1.1 ±0.1) X IQ- 


-3 


(1.7 ±0.3) X 10- 


-4 


(2.9 ± 1.5) X 10- 


-6 


2 ~ 2 


DE (o- = 0.36) 


(5.7 ±0.5) X IQ- 


-4 


(1.7 ±0.2) X 10- 


-4 


(2.0 ±0.4) X 10- 


-5 




Schechter 


(5.8 ±0.4) X 10- 


-4 


(1.7 ±0.2) X 10- 


-4 


(1.7 ±0.4) X 10- 


-5 



TABLE 6 

The REST-FRAME 8 fim LF for all 

GALAXIES AT 2 ~ 2 OBTAINED WITH THE 
1/Vmax METHOD. 



log^oi^Lt"'^) logio*(Mpc-^dcx-i) 



11.09 


o r,-, + 0.05 

■J-^°-0.09 


11.29 


_o 41 +0.08 


11.49 


n 00+0. 18 
■^■"■^-0.12 


11.69 


-4.28l°-| 


11.99 


-4.88^^3! 



TABLE 7 

Parameter values characterising the bolometric IR LF for star-forming galaxies 

AT 2 = 1 and 2 ~ 2. 

Redshift Functional form a a Ljj^(Lq) (Mpc^'^dex^-'^) 

2 = 1 Double-exp (eq. [TJ 1.2 0.39 (fixed) (2.5^0.3) x lO" {i-^tl'l) X 10"^ 
2 ~ 2 Double-exp (eq. [ij 1.2 0.39 (fixed) (6.3^0 9) X lO" (9.2^^ 7) X 10""' 



TABLE 8 

Number densities of star-forming LIRG and ULIRG at different redshifts. These 
number densities have been obtained by integrating the functional form 

APPEARING in THE SECOND COLUMN AND ARE EXPRESSED IN UNITS OF Mpc~^. DE STANDS 

FOR DOUBLE-EXPONENTIAL. 



Redshift Functional form logio(iboL) > H LIRG ULIRG 

2~0 DE(cr = 0.39) (4.1 ±0.3) X 10""* (4.1 ± 0.3) x 10"* (3.9 ± 0.7) x 10"^ 

2 = 1 DE(cr = 0.39) (2.6 ± 0.1) X 10"^ (2.5 ± 0.2) x 10"^ (1.2 ± 0.2) x lO"* 

2~2 DE(cr = 0.39) (1.1 ±0.1) X 10-3 (9.5 ± 1.5) x 10-* (1.5 ± 0.2) x 10"* 
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Fig. 1. — The comparison between photometric and spectroscopic redshifts for galaxies in our 24 /ini -selected sample in the GOODS 
fields. The distribution of relative errors dz = (zpi^ot ^ Zspec)/(1 + Zspec) has a median -0.007 and a dispersion ctz = 0.05. The horizontal 
lines separate the galaxies with 1.7 < Zpi^ot < 2.3 and the vertical lines, those with 1.7 < Zspec < 2.3. The distribution of relative errors 
for the 1.7 < Zgpec < 2.3 subsample of galaxies has a median -0.01 and a dispersion <t^ = 0.06. 
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Fig. 2. — The distribution of rest-frame 8 /im luminosities for galaxies at redshifts 0.9 < z < 1.1 (left-hand panel) and 1.7 < z < 2.3 
(right-hand panel). In each panel, the empty and shaded histograms include all the galaxies and only those with spectroscopic redshifts, 
respectively. 
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Fig. 3. — i?-band magnitudes versus soft X-ray fluxes for the X-ray detected galaxies in our 24 /xm galaxy sample in the GOODS/HDFN. 
The error bar for a generic source with soft X-ray flux 10"^'' crgcm~^s~^ and R = 24 mag is shown. The left-pointing arrows indicate 
that the soft X-ray flux 3 x 10"'^'^ ergcm~-^s~^ is an upper limit. 
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Fig. 1. IRAC-bascd colour-colour diagram for the 24 /im sources with rcdsliifts z > 1.5 in the GOODS fields. Filled squares and empty 
circles refer to X-ray classified AGN and to all the other z > 1.5 24 /im galaxies, respectively. The cross-like symbols indicate the few 
star-forming galaxies a.t z > 1.5 which are X-ray detected. The typical error bars for the colours of these sources are indicated in the lower 
right corner of the plot. 
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Fig. 5. — The k-corrcctions between 11.4-12.7 and 8 ^tm fluxes obtained using different IR galaxy model templates; Lagache et al. (2004; 
solid lines); Chary & Elbaz (2001; dotted lines) and Dale & Helou (2002; dashed line). The thin solid (dotted) lines indicate the interval 
of corrections obtained using the different models of Lagache et al. (Chary & Elbaa) with bolometric IR luminosity Ljr > lO'^'^ Lq. The 
corresponding thick lines indicate median k-corrections. The dashed line corresponds to the Dale & Helou model with parameter a = 1.1 
and 1.4. 
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Fig. 6. — Upper panel: The rest-frame 8 fim LF for star-forming galaxies at z = 1 in the GOODS fields, compared to the 8 fim LF at z ~ 0. 
The cross-like symbols show the 8 fimLF for star-forming galaxies at z ~ 0, as computed by Huang et al. (2006) with the 1/Vmax method. 
The dotted line represents the best X'^"fit obtained using a double-exponential function as that in eq. |(TJ. The upward-pointing triangles 
show the 1/Vmax LF at z = 1 obtained in this work, only strictly in the region of completeness of 8 /.tm luminosities. Lines of different 
styles show the 8 fimLF at z = 1 computed with the ML STY analysis, assuming different laws: a double-exponential form with bright-end 
slope fixed to the local value (a = 0.36; solid line); the same double-exponential form with a free cr parameter (dotted-dashed line), and a 
Schechter function (dashed line). Lower panel: the 68.3 and 95.4% confidence levels in (tr, uLl ^™) space in the case of a double-exponential 
law with (T as a free parameter. The parameter values yielding the ML are a = 0.20~^Q'Qy and L* = i/L^ = (1.32^J'i^^) X 10^^ Lq. 
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log^Q[v (8 )im)] 

Fig. 7. — The rest-frame 8 /imLF for star-forming galaxies at ^ ~ 2 in the GOODS fields. The filled circles show the LP in the region of 
completeness of 8 /^m luminosities, as computed with the 1/Vmax method. The solid and dashed lines show the 8 /imhF at z ~ 2 computed 
with the ML STY method, assuming a double-exponential form as in eq. JTJ and a Schechter function, respectively. The filled square is an 
extension of the LF at the faint-end, obtained using stacking analysis (see text for details). The addition of this point a posteriori allows 
to validate the extrapolated shape of the LF at the faint end. The 8 firahF at 2: ~ computed by Huang et al. (2006) has also been added 
for a comparison. 
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log^Q[v (8 |im)] 

Fig. 8. — The compared rest-frame 8/xmLF for star-forming galaxies at z = 1 and z ~ 2, both obtained in the GOODS fields. Symbols 
and lines are the same as in Figure[7] The upward-pointing triangles correspond to the LF at z = 1, as computed with the 1/Vmax method. 
The dotted-dashed line is the result of the ML analysis at the same redshift, adopting a double-exponential law with a = 0.36. 
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Fig. 9. — The rest-frame 8/imLF for star-forming galaxies in the GOODS/CDFS and HDFN separated, as computed with the l/V^nax 
method. Left-hand panel: z = 1; right-hand panel: 2 ~ 2. 
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log^Q[v (8 |im)] 

Fig. 10. — The rest-frame 8 fimLF for all the 24 /xm -selected galaxies (i.e. star-forming galaxies and AGN) at z ~ 2. The diamond-like 
symbols indicate the LF computed with thcl/Vmax method. The dashed and solid lines show the LF computed with the ML analysis, 
assuming the functional form given in eq. I|10| l with a = 1.2 and 1.4, respectively. The remaining symbols and line styles are the same as 
in Figure |8] 
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Fig. 11. — The bolometric IR- versus rest-frame 8 /um -luminosity relation for galaxies with i/L^'""> 10^" Lq in the Bavouzet et al. (2006) 
sample. The cross-like symbols indicate individual galaxies at redshifts 0.0 < z < 0.6. The dashed line shows the best-fit relation. The 

filled circle shows the resulting average value of (j/LJ^™; ) for a sample of galaxies at 1.3 < z < 2.3, as obtained through stacking 

analysis in the GOODS/CDFS. This point shows that the average relation between vL^^^and L]^i for 1.3 < z < 2.3 galaxies is basically 
the same as for galaxies at 0.0 < z < 0.6. 




Fig. 12. — Comparison between different Lj^^j versus i^L^ '""relations (left-hand panel) and derived conversion factors versus i^L^ (right- 
hand panel), as obtained from different calibrations available in the literature. 
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Fig. 13. — The evolution of the bolometric IR LF for star-forming galaxies from redshift z = to 2 ~ 2. Symbols and line styles are the 
same as in Figure [8] 
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Fig. 14. — The bolometric IR LF obtained in this work compared to the determinations of other authors at similar redshifts; 2; ~ 1 
(left-hand panel) and z ~ 2 (right-hand panel). 
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Fig. 15. — The evolution of the comoving bolometric IR luminosity density with redshift. The filled upward-pointing triangle and 
circle at redshifts 2 = 1 and z = 1.93 indicate the estimations of the respective bolometric IR luminosity density obtained in this work: 
QiR = (1.2 ± 0.2) X 10^ and (6.6^10) x IO^LqMpc-^. The density at z = 0.2 has been obtained from the bolometric IR LF derived 
from the 8/imLF by Huang et al. (2006). The red thick solid line corresponds to an interpolation between these redshifts, assuming a 
[(1 + 22)7(1 + 21)]^ evolution. The red thin solid lines indicate error bars on this evolution. Blue dashed and and orange dotted-dashed lines 
show the contributions of LIRG and ULIRG, respectively, at different redshifts. Other symbols refer to IR luminosity densities taken from 
the literature, and based on different datasets; ISO mid-IR (Flores et al. 1999; left-hand-pointing triangles), Spitzer mid-TR (Le Floc'h et 
al. 2005 and Perez-Gonzalez et al. 2005; right-hand-pointing triangle and asterisks, respectively), sub-millimetre (Barger et al. 2000 and 
Chapman et al. 2005; small and large diamond-like symbols, respectively) and radio (Haarsma et al. 2000; downward-pointing triangle). 
Some of these IR luminosity densities have been obtained from the star-formation rate densities compiled by Hopkins et al. (2004) and 
converted with the Kennicutt (1998) formula SFR = 1.72 X IQ-^" Lir. 



